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ABSTRACT 


The  eight-county  Allegheny  Region  of  Pennsylvania  Is 
84 5C  forested  and  comprises  one  of  the  state's  most  valuable 
tlmbersheds.  Much  of  this  area  Is  Inaccessible  to 
conventional  logging  equipment  due  to  steep  slopes,  poor 
drainage,  and  high  erodablllty.  For  this  reason,  there  has 
been  a  recently  growing  Interest  In  using  cable  yarders  to 
harvest  these  areas. 

A  computer  simulation  model  was  developed  to  aid  In 
Investigating  the  feasibility  of  cable  logging  on  a 
particular  site  In  the  region.  Data  used  were  those 
collected  from  a  cable  logging  operation  that  was  conducted 
on  the  Allegheny  National  Forest  In  October  of  1986. 

The  simulation  model  utilizes  the  SIMAN  simulation 
language  and  simulates  the  functional  elements  of  a  cable 
yarding  operation.  The  utility  of  the  model  was 
highlighted  by  comparing  two  plausable  harvest  unit 
configurations  for  a  given  site.  The  alternative  with  the 
shortest  total  harvesting  time  was  then  chosen. 

Several  design  features  will  aid  In  the  expansion  of 
the  model  as  more  detailed  data  are  collected.  In 
particular,  the  model  Is  modularized  by  harvesting 
function,  contains  complete  1 1 ne-by-1  1  ne  documentation  of 
the  SImAN  source  code,  and  contains  a  thorough  discussion 
of  the  statistical  methodologies  used.  Recommendations 
were  made  for  further  data  acquisition  and  model 


refinement. 
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INTRODUCTION 


Study 


The  eight-county  Allegheny  Region  (Figure  1.1) 


Is  located  along  the  northern  tier  of  Pennsylvania  and 


comprises  one  of  the  state's  most  valuable  tlmbersheds.  The 


region  Is  84*  forested  and  contain*  the  highest  per  acre 


and  total  growing  stock  volume  of  any  region  In 


Pennsylvania  (Powell  and  Consldlne  1982).  It  contains  3.3 


million  acres*  over  20*  of  the  state's  commercial  forest 


land  (Lord  1985).  Most  of  this  Is  of  the  valuable  northern 


hardwood  forest  type  consisting  of  high  quality  black 


cherry,  white  ash,  and  sugar  maple. 


Much  of  this  acreage  Is  characterized  by  steep  slopes. 


poor  drainage,  and  high  erodlblllty.  On  the  Allegheny 


National  Forest  alone.  It  Is  estimated  that  25,000  to 


40,000  acres  are  eligible  for  cable  yarding,  due  to 


drainage  and/or  slope  problems  (Hocklnson  1986). 


Environmental  concerns  over  high  soil  erosion  and  stream 


sedimentation  have  precluded  the  use  of  conventional 


logging  equipment,  such  as  rubber-tired  skldders,  to  remove 


timber  from  these  sites. 


Cable  yarding  systems  were  used  to  skid  timber  In  the 


eastern  United  States  In  the  early  1900's.  Their  use  died 


out  as  the  old  growth  stands  disappeared.  In  the  1970's, 


cable  yarders  were  reintroduced  to  harvest  steep  terrain 


m 


ALLEGHENY  REGION 


and  environmentally  sensitive  areas  (Peters  1984).  Cable 
logging  has  a  less  detrimental  Impact  on  the  environment 
(primarily  by  reducing  the  number  of  haul  and  skid  roads) 
than  conventional  systems.  For  this  reason  there  has  been 
a  recently  growing  Interest  In  using  cable  yarders  to  log 
environmentally  sensitive  areas  within  the  Allegheny 
National  Forest. 

Using  cable  logging  to  harvest  eastern  hardwood  logs 
on  steep  terrain,  however,  can  result  In  low  production 
rates  and  high  costs  per  unit  of  wood  produced  (LeDoux 
1985).  Logging  managers  can  Improve  productivity  and 
profitability  by  knowing  how  s 1 te-spec 1 f 1 c  variables 
Interact  with  cable  logging  equipment.  Carrying  out  this 
Investigation  by  field  study  or  trial  and  error  alone  can 
often  be  an  expensive  proposition.  One  tool  that  can  aid 
In  the  decision-making  process  Is  computer  simulation. 
Computer  simulation  provides  the  logging  manager  with  a  low 
cost  means  of  exploring  various  system  alternatives  before 
they  are  carried  out  and  thus  provides  a  valuable  tool  In 
the  decision-making  process. 

Ob  1 ect Ives 

This  study  Is  the  first  phase  In  the  development  of  a 
simulation  model  that  logging  planners  could  use  to 
Investigate  the  feasibility  of  cable  logging  on  a 
particular  site  and  under  particular  operating  conditions. 
The  basic  model  proposed  here  does  not  examine  the 
Interaction  between  logging  equipment  and  site  specific 


variables  but  can  be  refined  to  do  so  as  more  data  are 
collected  In  the  future. 

In  October  of  1986,  a  cable  logging  operation  was 
performed  on  the  Allegheny  National  Forest  to  Investigate 
the  feasibility  of  its  use  on  a  large  scale  In  the  region. 

The  data  used  In  the  development  of  this  model  are  limited 
to  those  obtained  from  this  one  time  and  motion  study.  The 
small  data  base  will  limit  the  actual  use  of  this 
particular  model  to  plan  cable-logging  operations  on  the 
Allegheny  National  Forest. 

The  specific  objectives  of  this  study  were  the  following 

1)  design  a  model  that  simulates  the  cable  yarding 
operation  that  was  performed; 

2)  construct  such  a  model,  using  SIMAN,  with  field  data 
that  was  collected  during  the  operation; 

3)  evaluate  the  model  with  regard  to  several  criteria. 
Including  acceptability  of  model  design,  appropriateness 
of  regression  equations,  theoretical  distributions,  and 
empirical  distributions  used,  and  ease  of  expanding  the 
model  as  more  data  are  collected; 

4)  make  recommendations  for  further  data  acquisition  and 
model  refinement. 


Chapter  II 


REVIEW  OF  LITERATURE 

Introduction 

This  chapter  will  present  a  review  of  the  literature 
on  the  following  topics:  (1)  how  elemental  cycle  times 
of  skyline  logging  operations  have  been  defined  In  the  past 
for  eastern  harvesting  operations,  (2)  the  development  of 
timber-harvest  simulation  models  In  the  past,  and  (3) 
the  validation  and  verification  of  simulation  models.  This 
chapter  will  serve  primarily  as  background  and  support  for 
the  methods  to  be  used  In  the  study. 

Elemental  Cvcl o  Times 

Since  1971,  the  USDA  Forest  Service  Engineering 
Research  Unit  In  Morgantown,  West  Virginia  has  been  active 
In  studying  eastern  cable  logging  operations.  As  a  result, 
turn  time  predicting  equations  for  various  cable  yarding 
systems  are  abundant  and  well  documented  In  a  number  of 
publications:  Cubbage  and  Gorse  (1975),  Fisher  et  al . 
(1980),  Rossle  (1983),  Biller  and  Peters  (1982  and  1984), 
Peters  (1984),  Fisher  et  al.  (1984),  Biller  and  Fisher 
(1984),  Peters  and  Baumgras  (1984),  Baumgras  and  Peters 
(1985),  LeDoux  (1985),  and  LeDoux  and  Starnes  (1986). 

In  all  of  the  studies,  total  turn  time  was  broken  down  Into 
five  elements,  and  delay-free  elemental  yarding  prediction 
equations  were  developed.  These  five  elements  were 
outhaul,  hooking,  lateral  Inhaul,  Inhaul,  and  unhooking. 
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Additionally,  a  delay-free  total  cycle  time  equation  was 
provided  In  each  of  the  studies.  This  equation  was 
developed  Independently  of  the  five  elemental  equations  and 


was  not  simply  the  sum  of  the  parts.  Since  some  of  the 
Independent  variables  were  highly  correlated,  they  would 
not  all  be  needed  to  predict  total  cycle  time.  Thus,  an 


Independent  equation  was  developed. 

Linear  Harvest  simul at 


Model s 


Many  different  1 1  mb e r- ha rvest  simulation  models  have 
been  developed  over  the  past  two  decades,  and  they 
represent  a  number  of  modeling  viewpoints.  Some,  such  as 
those  presented  by  Stark  (1968),  Bussell  et  al .  (1969), 
Amer Ican-Pul pwood  Association  (1972),  Martin  (1975), 
Klllham  (1975),  and  Johnson  (1976)  model  a  variety  of 
systems . 

Others  pertain  to  specific  systems:  Johnson  (1970) 
modeled  the  loading  and  hauling  subsystems  of  a  logging 
system;  Johnson  and  Biller  (1973)  modeled  a  wood  chipping 
system;  Bradley  et  al .  (1976)  and  Bradley  and  Wlnsauer 
(1976)  modeled  a  whole  tree  field  chipping  operation  and; 
Bare  et  al .  (1976)  modeled  a  logging  residue  handling 
system. 

Goulet  et  al .  (1980)  evaluated  five  of  these  models: 
Simulation  Applied  to  Logging  Systems  (SAPLOS)  (Johnson 
1976),  Timber  Harvest  and  Transport  Simulator  (THATS) 
(Martin  1975),  Full  Tree  Field  Chipping  and  Transport 
Simulator  (FTFC)  (Bradley  et  al.  1976  and  Bradley  and 
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Wlnsauer  1976),  Forest  Harvest  Simulation  Model  (FHSM) 


(Klllham  1975),  and  Harvest  System  Simulator  (HSS)  (American 
Pulpwood  Association  1972).  They  found  that  while  many 
user  Implementation  problems  exist,  the  models  still 
present  a  good  picture  of  state  of  the  art  In  timber 
harvesting  computer  simulation.  The  salient  features  of 
each  of  these  models  are  presented  In  Appendix  A;  Goulet's 
conclusions  are  presented  below. 

Implementation  ErflLlajLS 

Goulet  et  al .  (1980)  installed  the  five  models 
on  Auburn  University’s  IBM  S  370/158  computer  and  ran  them 
with  test  data.  In  summarizing  their  study,  Goulet 
concluded  the  following: 

1)  Each  model  operates  under  a  slightly  different  set  of 
rules  and  assumptions  according  to  the  philosophy  of  the 
model  builder.  Users  are  advised  to  proceed  with 
caution  when  choosing  and  using  a  model  and  to  be  aware 
of  the  assumptions  made  so  that  output  can  be  analyzed 
In  this  light. 

2)  The  models  are  not  easy  to  use,  and  In  general,  close 
coordination  between  a  computer  specialist  and  the  user 
will  be  necessary. 

3)  The  models  FHSM,  FTFC,  SAPIOS,  and  THATS  would  require 
extensive  design  and  reprograming  to  simulate  systems  not 
covered  In  the  basic  model. 

Despite  these  problems  the  authors  concluded  that  the 
models  represented  a  very  good  picture  of  state  of  the  art 
In  timber  harvesting  simulation.  Furthermore,  much 
learning  occurred  In  the  generation  of  these  models  and 
many  problems  were  uncovered  In  modeling  and  model 
1 mpl ementat Ion. 


They  stressed  that  continued  development  and 

refinement  of  tlmDer  harvesting  simulation  models  were 

needed  to  effectively  analyze  current  and  proposed 

harvesting  strategies.  They  suggested  that  the  results  of 

the  research  and  the  learning  derived  from  the  present 

models  be  incorporated  in  a  nee  model  that: 

a)  Is  des'gned  and  written  from  the  user's  point  of  view; 

b  )  faithfully  reproduces  the  harvesting  operations  to  be 
modeled; 

c)  maintains  a  level  of  detail  that  Is  uniform  across  all 
/unctions; 

d)  collects  model  statistics  to  estimate  both  the  mean  and 
the  variance  of  each  performance  variable; 

e)  de^‘nes  performance  variables  which  can  be  used  to  study 
the  b a  1 ance/ 1  mb  a  1 ance  of  the  system,  measure  the  complex 
interaction  of  personnel  and  machines,  and  effectively 
measure  marginal  and  total  costs; 

f!  allows  flexibility  for  tailoring  to  existing  systems  and 
for  the  creat'on  of  new  systems;  and 

g)  Is  usable  without  extensive  computer  training. 

&  Second  Ganerat 1  on  Model 

Work  py  Hines  at  al.  (1981)  established  further  design 
cr1ter’'a  for  a  second  generation  harvesting  simulation, 
particularly  mo d u 1  a r i z a t 1 o n  by  harvesting  function  within 
an  overall  s'mulatlon  framework  established  by  the  SLAM 
(pr'isker  1984)  simulation  language.  The  modular  design 
strategy  proposed  by  wines  et  al.  (1981)  envisioned 
separate  modules,  or  ;  u  ‘  1  a  ’  -i  q  blocks,  each  modeling  a 
unique  harvesting  function.  A  simulation  run  wou’d  Involve 
assembl'ng  the  appropriate  modules  and  executing  a  control 


program  that  would  pass  resources  between  modules  and 
collect  model  statistics. 

To  date,  three  modules  have  been  completed.  Webster 
et  al.  (1983)  reported  the  development  and  general  features 
of  a  f el  1 er/bu ncher  module  that  simulates  the  operation  of 
one  or  two  small,  skid-steer  machines.  A  highly  detailed 
description  of  the  f el  1 e r/ b u nc her  module  Is  given  by 
Padgett  (1982).  Hines  et  al.  (1983)  reported  the  general 
features  of  a  grapple  skidding  module  that  can  model  one  or 
two  grapple  skldders.  Liu  (1981)  provided  a  detailed 
description.  An  Input  data  pre-processor  module  has  also 
been  developed  (Rogers  1984)  to  assist  the  user  In  entering 
data  for  the  fel 1 er/buncher  modules.  Personal 
communication  with  Rummer  (1986)  suggested  that  further 
model  refinement  was  necessary  before  the  model  could  be 
used  for  planning  actual  logging  operations. 

Yal 1 d at  1  on  and  Verification 

One  of  the  most  Important  phases  In  the  development  of 
a  computer  simulation  model  Is  determining  whether  the 
model  Is  an  accurate  representation  of  the  real-world 
system  being  studied.  Model  developers  address  this 
concern  through  model  verification  and  validation.  Law  and 
Kelton  (1982)  described  verification  as  the  determination  of 
whether  a  simulation  model  performs  as  Intended,  i.e., 
debugging  the  computer  program.  They  defined  validation  as 
the  determination  of  whether  a  simulation  model  Is  an 
accurate  representation  of  the  real -wor ' 0  system. 


A  model  should  be  developed  for  a  specific  purpose  or 
use,  and  Its  validity  should  be  determined  with  respect  to 
that  purpose.  A  model  may  be  valid  for  one  set  of 
experimental  conditions  and  Invalid  for  another.  A  model 
Is  considered  valid  for  a  set  of  experimental  conditions  if 
Its  accuracy  Is  within  the  acceptable  range  required  for 
the  model's  Intended  purpose  (Sargent  1984). 

Below  Is  a  description  of  some  of  the  validation 
techniques  (and  tests)  used  In  model  validation. 

1)  Face  validity:  Face  validity  Is  asking  people 
knowledgeable  about  the  system  whether  the  model  and/or 
Its  behavior  Is  reasonable.  This  technique  can  be  used 
In  determining  If  the  logic  In  the  model  flowchart  Is 
correct  and.  If  a  model's  Input-output  relationships  are 
reasonable  (Sargent  1984). 

2)  Comparison  to  other  models:  Computer  output  from  the 
simulation  model  being  validated  Is  compared  to  the 
output  of  other  (valid)  models  (Shannon  1975). 

3)  Turing  Tests:  People  who  are  knowlegeable  about  the 
operations  of  a  system  are  asked  if  they  can 
discriminate  between  real-world  system  data  and  model 
output.  A  statistical  procedure  for  Turing  Tests  Is  given 
1 n  Schruben  (  1980  )  . 

4)  Historical  methods:  Naylor  and  Finger  (1967)  proposed 
the  three  historical  methods  of  validation:  rationalism, 
empiricism,  and  positive  economics.  Rationalism  assumes 
that  everyone  knows  whether  the  underlying  assumptions  of 
a  model  are  true.  Then  logic  deductions  are  used  from 
these  assumptions  to  develop  the  correct  (valid)  model. 
Empiricism  requires  every  assumption  and  outcome  to  be 
empirically  validated.  Positive  economics  requires  only 
that  the  model  be  able  to  predict  the  future  and  is  not 
concerned  with  Its  assumptions  or  structure. 

5)  Traces:  The  behavior  of  different  types  of  spec  '  f  1  c 
entitles  In  the  model  are  traced  through  the  model  to 
determine  If  the  model's  logic  Is  correct  and  If  the 

ecesary  accuracy  Is  obtained. 
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Some  of  the  comments  on  validation,  found  In  the 

documentation  manuals  of  the  models  that  were  reviewed,  are 

Instructive  In  showing  the  approach  to  validation  that  has 

Oeen  taken  In  these  models,  and  so  will  be  repeated  here. 

Johnson  et  al.  (1972,  p.  361)  say  of  SAPLOS, 

The  model  needs  further  validation  from 
studies  of  a  variety  of  logging  systems 
In  different  locations.  This  validation 
process  will  be  performed  In  conjunction 
with  the  expansion  of  the  data  base. 

Martin  (1975  ,  p  .  3  1  )  says  of  THATS, 

The  methodology  of  the  model  and  its 
structure  need  little  validation  because 
the  simulator  simply  duplicates  the 
process  of  logging  as  It  Is  normally 
performed  In  Appalachia.  The  sequence 
of  operations  Is  the  same  In  both. 

Bradley  et  al.  (1976,  p.ll)  caution  the  user  of  FTFC, 

The  user  must  test  the  simulator  on  an 
actual  logging  operation.  The  test  is 
required  to  detect  errors  in  either 
concept  or  model. 

Martin's  (1975)  comment  suggests  that  he  used  face 
validity  In  the  validation  of  THATS.  from  the  1 Iterature 
reviewed.  It  Is  uncertain  If  any  formal  types  of  va 1  1  flat ( on 
were  used  In  the  development  of  the  two  other  models. 

Webster  et  a  1 .  (1984)  suggested  that  t  r  v  '  n  ;  tc  level  :p 

a  valid  model  Involves  a  multitude  c  f  compromises  w  * • c h 
embrace  design,  practical,  logical,  econom1  ca'  ,  c  o  m  p  u  t“r  , 


anguage.  Implementation,  and  phi  1 osopn  ’  ca  ’  d  '  ‘ 


t  '  e  s 


At  the  heart  of  the  d  rob  1 em  of  d  e  v  e  ’  o  p  {  n  q  s  '  m  u  ‘  a  t  '  o  n  mode' 
Is  the  question  of  complexity  of  the  mode'.  To  he 
realistic,  a  model  may  need  to  be  romp’o.,  s'  nee  by 


\  V  *  \ 
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Including  more  of  the  factors  and  variations  of  the  system 
the  model  becomes  more  capable  of  duplicating  the  system's 
response.  But  complex  models  take  more  time  to  develop* 
are  more  difficult  to  analyze*  require  more  computer 
resources*  and  cost  more  overall  than  simple  models. 

Simple  models*  on  the  other  hand*  will  not  always  provide 
the  user  with  the  quality  or  variability  of  results  that 
are  naturally  Inherent  within  the  system  and  may  lead  him 
Into  erroneous  conclusions. 


? 


✓ 

t 
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Verification 

Law  and  Kelton  (1982)  describe  five  techniques  which 

can  be  used  for  debugging  the  computer  code  of  a  simulation 

model  . 

1)  In  developing  the  simulation  model*  write  and  debug  the 
computer  program  In  modules  or  subprograms. 

2)  Structured-walk-through:  All  people  Involved  In  the 
model  development  are  assembled  In  a  room.  They  go 
through  the  computer  code  and  do  not  proceed  from  one 
statement  to  another  until  everyone  agrees  that  the 
statement  Is  correct. 

3)  Trace:  In  a  trace,  the  state  of  the  simulated  system, 
1.e.»  the  contents  of  the  event  list,  the  state 
variables,  certain  statistical  counters,  etc.,  Is 
printed  Just  after  each  event  occurs  In  order  to  see 
whether  the  program  is  operating  as  Intended.  These 
traces  should  be  examined  In  order  to  see  if  the  model 
logic  Is  correct.  Most  major  simulation  language 
packages  provide  the  capability  to  perform  a  trace. 

4)  The  model,  when  possible,  should  be  run  under 
simplifying  assumptions  for  which  the  expected  model 
output  can  easily  be  computed.  The  expected  output 
should  then  be  compared  to  the  actual  output. 

5!  with  some  types  of  simulation  models,  It  may  be  helpful 
to  display  the  simulation  output  on  a  graphics  terminal 
as  the  simulation  actually  progresses. 


f. 


-*a.  A/» 
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Chapter  III 


DESCRIPTION  OF  THE  SIMULATION  STUDY 
Introduction 

A  computer  model  to  simulate  the  harvesting  actions  of 
a  cable  yarder  In  the  Allegheny  Region  of  Northwestern 
Pennsylvania  was  constructed  using  field  data  collected 
from  that  region.  The  model  predicts  time  to  harvest  a 
given  site  and  volume  harvested  per  unit  time.  The 
discrete  event  model  employs  the  SIMAN  simulation  language. 
This  chapter  describes  the  "real-world"  system  that  was 
observed  and  the  simulation  model  that  was  constructed. 

Iilfi  "Real-World"  System 

The  focus  of  this  section  Is  the  description  of  the 
logging  site,  the  logging  machine,  the  logging  operation 
that  was  conducted,  and  the  time  study  that  was 
performed . 

IM  Ig-5.1  liLSA 

The  logging  operation  was  conducted  on  the  Bradford 
Ranger  District  of  the  Allegheny  National  Forest  In 
Northwestern  Pennsylvania.  The  total  area  harvested  was  II 
acres.  The  primary  tree  species  were  black  cherry,  white 
ash,  sugar  maple,  and  basswood.  The  clearcut  harvest 
yielded  an  average  of  12.3  MBF/acre  of  sawtlmber  and  9.7 
cords/acre  of  pulpwood  (Table  3.1),  contained  In  an  average 
of  66  trees/acre  (Table  3.2). 
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The  cable  machine  used  was  a  slacking  (live)  skyline 
(Figure  3.1)  with  a  Christy  carriage.  The  yarder,  used  to 
transport  whole  trees  from  the  stump  to  a  landing,  was  not 
a  commercial  design,  and  would  cost  about  $50,000  If  built 
to  order.  It  Is  owned  and  operated  by  Bess  Skyline  Logging 
of  Virginia.  It  was  a  single-drum  yarder,  powered  by  a  six 
cylinder  Industrial  Waukesha  gasoline  engine.  Tower  height 
was  40  feet;  the  skyline  was  3/4-1nch  IWRl  regular-lay  wire 
rope;  the  mainline  was  l/2-1nch  IWRL  regular-lay  wire  rope; 
maximum  mainline  pull  was  15,000  pounds. 

The  Yarding  Qperat 1  on 

The  layout  of  the  harvesting  area  Is  Illustrated  In 
Figure  3.2.  Two  log  landings  were  constructed  for  the 
operation  (see  Figure  3.2).  The  yarder  was  positioned  at 
the  first  landing.  The  tower  was  raised  and  guyed  to 
stumps  or  trees  above  the  landing.  The  skyline  was  tied 
off  to  a  talltree  downhill  from  the  yarder.  The  area 
from  which  logs  were  harvested,  at  any  one  of  these  particular 
setups.  Is  termed  a  logging  corridor.  Logs  were  harvested 
from  several  corridors  at  each  of  the  two  landings. 

Changing  logging  corridors  consisted  of  anchoring  the  skyline 
to  another  talltree  so  that  the  skyline  would  be  situated 
over  an  area  that  had  not  yet  been  harvested.  This 


particular  type  of  harvesting  pattern  Is  described  by 
Peters  (1985)  as  a  reverse  fan  type;  the  cable  yarder  being 


GUYLINE 


Schematic  diagram  of  a  cable  logging  operation. 


Logging  Corridors 


the  point  of  common  pivot  and  the  talltrees  creating  the 
arc  of  the  fan. 

After  a  corridor  was  set  up,  the  Christy 
carriage  traveled  down  the  tight  skyline  until  It  engaged  a 
stop,  which  released  the  mainline  from  the  carriage.  Three 
chokesetters  took  the  mainline  out  and  attached  the  logs, 
then  signaled  the  yarder  operator  to  winch  In  the  load. 

When  the  load  reached  the  carriage,  a  ball  on  the  mainline 
unlocked  the  carriage  from  the  stop.  The  load  traveled 
uphill  to  the  yarder  at  the  log  landing,  where  the  skyline 
was  slackened  and  the  logs  were  unhooked  by  a  single 
unhooker.  These  six  elements  ( outhaul , 1 ateral  outhaul, 
hooking,  lateral  Inhaul,  Inhaul,  and  unhooking)  comprise 
one  yarding  cycle.  For  each  corridor  change,  yarding 
started  at  the  top  of  the  hill  and  progressed  toward  the 
bottom.  Maximum  slope  yarding  distance  was  425  feet.  One 
stem  was  usually  landed  per  turn  and  occasionally  two  If 
the  trees  were  small. 

Logs  that  accumulated  at  the  landing  area  were 
Intermittently  transported  (by  means  of  a  rubber-tired 
skldder)  to  a  main  loading  area.  It  was  necessary  for  the 
cable  yarder  to  cease  operation  while  the  skldder  was  at 
the  landing. 

When  all  corridors  were  harvested  at  landing  1,  the 
guy  lines  were  disconnected,  the  tower  was  lowered,  and 
the  mainline  and  skyline  were  wound  on  the  drum.  The  cable 


yarder  was  then  transported  to  landing  2  where  It  was 
assembled  by  reversing  this  process. 

The  conditions  at  landing  1  differed  from  those  at 
landing  2  In  several  ways.  Landing  1  was  situated  In 
such  a  way  that  the  skyline  crossed  over  the  main  logging 
road  (see  Figure  3.2).  On  several  occasions*  yarding  was 
suspended  In  order  to  allow  a  logging  truck  to  pass;  this 
was  not  the  case  for  landing  2.  The  slopes  were  much 
more  gentle  at  landing  2.  This  caused  the  logs  to  drag  on 
the  ground  on  the  Inhaul  phase.  This*  combined  with  the 
fact  that  the  part  of  the  unit  harvested  from  landing  2 
was  much  brushler,  caused  logs  to  get  hung  up  In  the  slash 
more  often  on  the  Inhaul  phase.  Separate  service  time 
equations*  for  some  of  the  yarding  cycle  elements*  were 
developed  for  each  landing  In  order  to  explain  some  of  this 
var 1 ab 1 1 1 ty . 

Timing  the  Operation 

Time  and  motion  study  data  were  collected  over  a  two 
week  period  In  October  of  1986.  Statistics  were  collected 
on  183  yarding  cycles  from  five  corridors  at  two  landings; 
three  corridors  at  the  first  landing  and  two  at  the  second. 
Throughout  the  study,  continuous  timing  was  used  to 
document  elemental  cycle  times  and  delay  times  to  the 
nearest  one  hundredth  of  a  second. 

Most  of  the  sampling  was  devoted  to  timing  the 
productive  yarding  elements  so  that  elemental  time 


predicting  equations  could  be  developed  for  each  of  the  six 


phases  of  a  yarding  cycle,  or  turn  (outhaul,  lateral 
outhaul,  hooking,  lateral  Inhaul,  Inhaul,  and  unhooking). 
Additionally,  all  hauling  distances  as  well  as  number  of 
stems  yarded  per  turn  were  recorded. 

Four  non-productive  delay  times  were  also  recorded: 
time  to  change  landings,  time  to  change  corridors,  time 
needed  to  clear  the  landing  of  accumulated  logs,  and  time 
required  to  free  logs  hung  up  In  the  slash  on  the  Inhaul 
element.  The  statistics  collected  are  provided  In  Appendix 
C,  and  their  definitions  are  In  Table  3.3. 

IM  Slmul  at  1  on  Mg^gJ 

This  discussion  will  address  four  topics.  First,  the 
functional  elements  of  the  "real-world"  system  the  model 
considers  and  their  dynamic  relationship  will  be  presented. 
Two  flowcharts  will  be  employed  for  this  purpose:  the  Main 
flowchart  (Figure  3.3)  and  the  Yarding  Cycle  flowchart 
(Figure  3.4).  Attention  Is  then  directed  to  the 
development  of  random  Inputs  which  represent  service  times 
and  Incremental  distance  changes  In  the  model.  Following 
this.  Important  features  of  SIMAN  and  of  model  construction 
will  be  pointed  out.  A  section  on  model  validation  will 
follow. 

A  complete  executable  example  that  compares  two 
system  alternatives,  output,  and  directions  for  modifying 
the  model  (e.g.  such  as  changing  the  number  of  corridors, 
corridor  lengths,  etc.)  Is  provided  In  Appendix  B.  This 
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Table  3.3:  Definitions  of  model  statistics  collected. 

CPthaul  time  and  distance:  CLrthaul  ends  when  the  carriage  hits  the  stop. 

Lateral  outhaul  tine  and  distance:  LartBral  outhaul  ends  when  the  chokers  reach 
the  stars. 

Lateral  Inhaul  time:  Lateral  Inhaul  tins  ends  when  the  carriage  begins  to  move 
uphill. 

Hooking  tins:  Hocking  time  ends  when  the  signal  Is  give  to  haul  In. 

Inhaul  tlffB:  Inhaul  tine  ends  trtien  the  stems  are  on  the  ground  at  the  landing. 

Unhooking  time:  Unhooking  tins  ends  wnen  the  outhaul  begins. 

Total  cycle  tins:  Total  cycle  time  begins  and  ends  at  the  start  of  the  outhaul . 

Nurber  of  stars  landed:  This  Is  the  muter  of  stars  per  turn  that  are 
successfully  yarded  to  the  landing. 

Svstan  delays  Cfrecuencv  and  duration):  Systan  delays  are  defined  as  any  event 
that  disrupts  the  "normal  flow*’  of  activity.  Systan  delays  that  were  recorded 
are  as  follows: 

*  Delay  to  clear  landing:  This  tin©  begins  when  unhooking  ends  and  ends  when 
outhaul  begins.  The  duration  of  this  delay  is  0.00  unless  logs  are  actually 
cleared  from  the  landing  on  that  cycle. 

*  Delay  ic  change  corridors:  This  tire  begins  when  unhooking  ends  at  the 
previous  corridor  and  ends  when  outhaul  begins  at  the  new  corridor. 

*  Delay  1c  change  landings:  This  tire  begins  when  unhooking  ends  at  the 
previous  landing  and  ends  when  outhaul  begins  at  the  new  corridor. 

*  Logs  hung  in  slash  an  inhaul  element:  This  time  beings  when  carriage 

ret  Ion  stops  on  the  Inhaul  element  and  ends  when  motion  of  the  carriage  resumes 
H.e.  after  log  Is  freed  from  slash!. 


Figure  3.3:  Main  flowchart 


*  YARDING  CYCLE  FLOWCHART  ' 


Go  to 


on  mam  flowcrjit 


c'qwe  3  4:  Yarding  cycle  flowchart. 


example  serves  to  Illustrate  the  potential  usefulness  of  a 
model  of  this  type, 
lilfi  Malfl  Flowchart 

The  Main  flowchart  (Figure  3.3)  contains  the  model 
logic  that  represents  movement  between  corridors,  and 
between  landings.  The  static  and  dynamic  relationships  are 
1 1 1 ustrated . 

The  number  of  landings  In  the  harvesting  unit  are 
read  Into  the  model  and  a  delay  to  set  up  a  landing  Is 
executed.  The  number  of  corridors  at  the  current  landing, 
the  current  corridor  length,  and  Initial  outhaul  distance 
(all  specified  by  the  user)  are  set.  A  yarding  cycle  Is 
then  completed. 

The  number  of  cycles  to  be  completed  at  the  current 
outhaul  distance  Is  then  sampled  from  the  user-defined 
empirical  distribution.  When  all  cycles  are  completed  at 
that  distance,  outhaul  distance  Is  Increased  by  an 
Increment  that  was  specified  by  the  user. 

If  the  outhaul  distance  Is  less  than  the  current 
corridor  length,  more  yarding  cycles  are  performed.  If 
this  Is  not  the  case,  a  check  is  made  to  determine  If  there 
are  any  unharvested  corridors  at  the  current  landing.  If 
there  are,  a  delay  to  set  up  a  new  corridor  Is  executed  and 
the  corridor  Is  harvested  In  the  same  manner  as  described 
above . 

If  there  are  no  more  unharvested  corridors  at  the 


current  landing,  than  a  check  Is  made  to  determine  if  there 


are  any  more  landings  In  the  harvesting  unit.  It  there 
are,  a  delay  to  set  up  the  new  landing  is  executed;  if  not, 
the  simulation  Is  terminated. 


ItLfl  La r a i n a  Cxsla  Elflj&HAit 

The  Yarding  Cycle  flowchart  (Figure  3.4)  contains  the 
model  logic  that  represents  a  yarding  cycle.  The  static 
and  dynamic  relationships  are  Illustrated. 

The  six  productive  elements  of  a  caple  yarding  cycle 
have  been  considered;  outhaul,  lateral  outhaul,  hooking, 
lateral  Inhaul,  Inhaul,  and  unhooking.  Additionally,  two 
nonproductive  delays  have  been  considered:  Inhaul  delays, 
and  a  delay  time  to  clear  the  landing  of  accumulated  logs. 
Two  other  nonproductive  delays  (delay  to  change  corridors 
and  to  change  landings)  are  Included  in  the  Main  flowchart. 

ca.l  Aralys f.s 

This  section  describes  how  the  data  collected  were 
used  to  specify  random  Inputs  In  the  simulation  model. 
Theoretical  distributions,  empirical  distributions, 
regression  equations,  and  median  values  were  used. 

Thsgrat  1  cal  D  1  str  i  but  i  ons 

Standard  techniques  of  statistical  Inference  were  used 
to  fit  a  theoretical  distribution  to  the  data  for  some  of 
the  cycle  elements.  After  a  theoretical  distribution  was 
hypothesized,  maximum  likelihood  estimates  of  the 
distribution  parameters  were  calculated  using  the 
Statgraphics  statistical  package  (STSC,  Inc.  1985).  The 
Chi-sauare  goodness-of-f it  test  and  K o 1  mo ao ro v - S m 1 r n o v  test 


•ere  employed  to  determine  If  an  acceptable  fit  was 
obtained.  Random  numbers  were  sampled  directly  from  these 
theoretical  distributions  In  the  simulation.  Probability 
density  functions,  of  the  theoretical  distributions  used, 
are  provided  In  Table  3.4. 

Service  time  distribution  results  for  hooklpg  (Figure 
3.5  and  Table  3.5),  unhooking  (Figure  3.6  and  Table  3.6), 
lateral  outhaul  for  landing  1  (Figure  3.7  and  Table  3.7), 
and  lateral  outhaul  for  landing  2  (Figure  3.8  and  Table 
3.8)  are  provided. 
tnuLirliil  BlxLclMilaiis 


When  a  theoretical  distribution  form  could  not  be 
found  to  adequately  fit  the  data,  the  data  were  used 
directly  to  define  an  empirical  distribution.  In  the 
simulation,  random  numbers  were  sampled  directly  from  this 
empirical  distribution.  This  method  was  employed  to 
describe  lateral  Inhaul  time  (Table  3.9)  and  number  of 
yarding  cycles  to  be  performed  at  a  given  outhaul  distance 
(Table  3.10). 

Reg  ress 1  on  Equations 

When  two  yarding  cycle  elements  were  highly 
correlated,  their  relationship  was  described  through 
standard  regression  techniques.  A  linear  relationship  was 
found  to  exist  between  outhaul  distance  and  outhaul  and 
Inhaul  service  times.  In  the  simulation,  outhaul  distance 
Is  used  to  predict  outhaul  time  and  Inhaul  time  In  two 
separate  regression  equations.  Outhaul  distance  Is  highly 


laDle  J.4:  KrooaDinty  aensity  runctions  ot  two  tneore 
distributions. 


Wei  bull  Distribution 

f(x)  •  afi"a  x3'1  e-(x/B)“  x>0 


Lognormal  Distribution 


f(*) 


e-(ln  x  - u)2/2cr2 
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HOOKING  TIME  (SECONDS) 

Figure  3.5:  Frequency  distribution  of  hooking  time  with 
fitted  Wei  bull  curve  superimposed. 
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Figure  3.7:  Frequency  distribution  of  lateral  outhaul  t 
for  landing  1  with  fitted  Weibull  curve  superimposed. 
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for  landing  2  with  fitted  Weibull  curve  superimposed. 
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correlated  to  Inhaul  time  because  the  value  of  outhaul 
distance  will  always  be  equal  to  Inhaul  distance. 

Regression  results  for  Inhaul  time  for  landing  1  and 
landing  2  (Table  3.11  and  Table  3.12  respectively)  and 
outhaul  time  for  landing  1  and  landing  2  (Table  3.13  and 
Table  3.14  respectively)  are  provided. 

Non-oarametr 1 c  Methods 

When  there  was  a  low  sample  size  on  a  particular 
element  of  Interest,  these  elements  were  represented  In  the 
simulation  model  by  using  the  median  value  of  the  data 
points  that  were  collected.  Median  values  were  used  to 
explain  the  frequency  and  duration  of  the  four 
nonproductive  delays  that  were  modeled.  Results  for  Inhaul 
delays,  delays  to  clear  the  landing,  delays  to  change 
corridors  and  to  change  landings  are  provided  In  Table 
3.15. 


Sel  ect Ion  al  1M  Language 

The  following  criteria  were  used  In  the  selection  of  a 
simulation  language: 

1)  The  completed  model  must  run  on  a  microcomputer. 

2)  The  language  must  possess  c ha racte r 1 st 1 cs  that  will  aid 
In  model  modification. 

3)  The  language  must  allow  the  collection  of  statistics 
during  the  simulation  run. 

4)  The  language  must  be  flexible  enough  to  adequately  model 
the  system  under  study. 

The  SIMAN  simulation  language  (Pegden  1985)  meets  all 
these  criteria.  It  Is  a  combined  discrete-continuous  event 


simulation  analysis  language  for  modeling  general  systems. 
Developed  In  1982  and  under  constant  revision,  it  Is 
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considered  to  be  state-of-the-art  In  current  simulation 


language  technology  (Klelndorfer  1986).  To  date  no  harvest 
simulators  exist  that  utilize  this  simulation  language. 
Discussion  here  will  be  limited  to  the  discrete  modeling 
capabilities  of  SIMAN. 

Structu  re 

SIMAN  Is  designed  around  a  logical  modeling  framework 
In  which  the  simulation  program  is  decomposed  Into  a  model 
frame  and  an  experimental  frame.  External  to  this,  SIMAN 
has  an  OUTPUT  processor  which  collects  and  analyzes 
simulation  results.  Debugging  the  model  code  Is  aided  by  a 
system  trace  and  Interactive  debugger. 

Model  F  rame .  The  model  frame  defines  the  static 
and  dynamic  characteristics  of  the  system.  Within  the 
model  frame,  either  an  event  or  process  orientation  can  be 
used  to  describe  the  model.  The  primary  modeling 
orientation  for  discrete  change  systems  Is  the  process 
orientation.  In  which  the  model  Is  constructed  by  depicting 
the  functional  operations  of  the  system  as  block  diagrams. 
The  block  diagram  Is  a  linear  top-down  sequence  of  blocks 
which  represents  specific  process  functions  such  as  time 
delays  and  queues.  This  orientation  was  used  In  the 
development  of  the  cable  logging  model. 

A  second  modeling  orientation  is  the  event  orientation 
which  may  be  used  to  augment  or  replace  the  block  diagram 
component  of  the  model  frame.  The  event  component 
consists  of  a  set  of  user  written  FORTRAN  subroutines  which 


contain  the  mathematical  and  logical  expressions  that 
define  Instantaneous  state  transitions  within  the  system. 

In  future  model  refinement  the  event  orientation  will 
be  used  to  augment  the  model  that  has  been  developed. 

Block  functions  (such  as  delay  for  hooking)  will  be  replaced 
perhaps  by  a  very  complex  subroutine,  thus  increasing  the 
level  of  detail  modeled. 

In  a  sense,  each  block  which  explains  a  delay  element 
of  the  logging  cycle  may  be  thought  of  as  a  module  In  the 
model.  Each  module  may  be  refined  Individually  by  simply 
replacing  a  specific  block  with  a  FORTRAN  subroutine  and 
plugging  It  back  Into  the  model. 

Exper Imental  Frame.  The  SIMAN  experimental  frame 
defines  the  experimental  conditions  under  which  the  model 
Is  to  be  run  In  order  to  generate  specific  output  data. 

This  Includes  such  elements  as  the  Initial  conditions  for 
the  run,  machine  capacities,  the  type  of  statistics  to  be 
recorded,  and  the  various  parameters  and  coefficients  of 
the  theoretical  d  1  st r 1 b u 1 1  on s ' and  regression  equations  that 
have  been  developed.  Since  these  elements  are  specified 
external  to  the  model  description,  they  may  be  easily 
changed  without  affecting  the  basic  model  definition.  Many 
different  scenarios  may  be  modeled  by  modifying  the 
experimental  frame. 

Output  Processor.  Based  on  the  model  and  experiment. 


the  SIMAN  simulation  program  generates  output  files  which 
record  the  model  state  transitions  as  they  occur  in 


simulated  time.  The  data  In  the  output  files  can  then  be 
subjected  to  various  data  analysis  within  the  SIMAN  output 
processor  or  exported  to  a  statistical  package  for 
analysis.  Within  the  SIMAN  framework,  the  data  analysis 
follow  the  development  and  running  of  the  simulation 
program  and  are  completely  distinct  from  It. 

System  Debugging.  Dobugging  or  model  verification  Is 
the  process  of  Isolating  and  correcting  the  logic  errors 
that  produce  Invalid  results.  The  SIMAN  system  trace  Is 
used  within  a  discrete  model  to  generate  a  detailed  trace 
report  of  the  processing  of  entitles.  In  the  event  mode, 
the  trace  report  summarizes  the  occurrence  of  each  event 
and  details  all  operations  executed  within  the  event.  When 
a  logic  error  Is  detected,  the  SIMAN  interactive  debugger 
may  be  used.  It  allows  the  user  to  Interactively  monitor 
and  control  the  execution  of  a  simulation.  Errors  can  be 
Isolated  and  corrected  during  execution  without  the  need  to 
recompile,  relink  and  rerun  the  simulation. 

Model  Val 1 d  at  1  on 

The  topic  of  validating  simulation  models  has  been 
discussed  In  the  validation  section  of  Chapter  II  and  so 
shall  not  be  repeated  here.  It  is  the  Intent  of  this 
discussion  to  describe  two  methods  used  In  the  validation 
of  this  model. 

It  was  discussed  earlier  that  a  model  should  be 
developed  with  respect  to  a  specific  purpose  and  validated 


for  that  purpose.  From  prior  discussion  In  this  chapter, 


should  be  apparent  that  the  purpose  of  this  model  was  to 
construct  a  working  model  that  will  pave  the  way  for  the 
development  of  a  more  detailed  model  in  the  future;  to  be 
used  as  a  guide  to  aid  further  simulation  studies.  The 
model  was  validated  with  respect  to  that  purpose. 

First,  the  model  was  examined  to  determine  If  It  had  a 
high  degree  of  face  validity.  The  model  flowcharts  (Figure 
3.3  and  Figure  3.4)  were  examined  to  determine  if  the  model 
logic  "mirrors"  the  logic  of  the  system  under  study.  Any 
simulation  study  embraces  a  series  of  compromises  as  to  the 
level  of  detail  that  should  be  modeled.  The  model 
developed,  contains  all  of  the  major  functional  elements  of 
the  system  studied  and  the  static  and  dynamic  relationships 
seem  to  make  sense.  Therefore,  the  model  arguably  contains 
a  high  degree  of  face  validity. 

The  following  procedure  was  used  to  validate  the  model 
empirically.  Each  of  the  theoretical  Input  distributions 
were  examined  for  goodness-of-f  i  t .  The  Chi-square  and 
Kolmogorov  Smirnov  goo d ne s s-o f- f  1 1  tests  were  used. 
Regression  equations  were  examined  and  detailed  residual 
analysis  was  performed.  Empirical  distributions  and 
median  values  are  simply  representative  of  the  data  points 
sampled.  A  detailed  discussion  of  these  statistical  results 
were  provided  In  the  statistical  analysis  section  of  this 
chapter  and  will  not  be  repeated  here. 


Chapter  IV 


SUMMARY  AND  CONCLUSIONS 

Sundry 

This  project  has  resulted  In  the  development  of  a 
timber  harvest  simulator  that  models  the  actions  of  a  cable 
yarder  operating  In  the  Allegheny  Region  of  Pennsylvania. 
Elemental  yarding  time  predicting  equations  were  developed 
from  field  data  collected  from  that  region. 

In  the  model  presented  in  this  study,  the  level  of 
detail  Is  the  function  itself  (e.g.  hooking  service  time, 
unhooking  service  time,  etc.).  Indirectly,  many  different 
"what-1f"  questions  may  be  answered,  but  It  requires  the 
user  to  estimate  the  change  In  the  maximum  likelihood 
estimates  of  the  parameters  affected  by  the  change 
proposed.  This  Is  Illustrated  In  Appendix  B  where  the 
model  Is  used  to  examine  two  harvest  unit  configurations 
In  order  to  choose  the  best  alternative  with  minimization 
of  make-span  as  the  performance  criteria. 

The  simulation  model  has  been  designed  around  the 
central  premise  that  a  simulation  model  will  go  through 
many  stages  of  development  during  its  lifecycle.  At  each 
stage,  the  model  will  be  refined  as  more  data  are 
collected;  perhaps  by  someone  that  has  not  been  involved  in 
the  study  to  date. 

For  this  reason  the  model  has  been  fully  documented. 
Each  line  of  SIMAN  code  has  been  commented  so  that  the 


model  logic  may  easily  be  understood.  The  methodology,  anc 
data  used  In  the  development  of  Input  distributions  have 
been  Included  and  the  results  Illustrated. 


As  previously  pointed  out,  the  model  Is  modular  In 
design.  Specific  functions  of  the  model  may  be  developed 
Independently  of  others.  Predicting  equations  developed 
from  the  time  study  conducted  may  be  used  for  functions 
where  further  data  are  not  available.  As  further  data  Is 
collected  on  other  functions,  these  functions  may  be 
developed  In  greater  detail.  Independently  of  the  others. 

Areas  Far.  Further  St.u.fl^ 

As  previously  mentioned,  at  the  heart  of  any 
simulation  study  Is  the  question  of  the  level  of  detail  to 
be  modeled.  Very  complex  models  take  more  time  and 
resources  to  develop  and  are  more  expensive  to  use.  Less 
complex  models  on  the  other  hand,  may  not  be  detailed 
enough  to  answer  specific  questions  of  Interest. 

The  approach  that  has  been  taken  In  the  model 
developed  In  this  study  was  to  model  the  elemental  yarding 
cycle  function  as  the  lowest  level  cf  detail.  Additional 
data  are  needed  If  a  higher  degree  of  detail  Is  to  Se 
modeled.  These  data  could  be  Incorporated  Into  a  new  model 
that  would  consider  specific  elements  w'thln  each  function 
that  has  been  currently  modeled.  Severa1  suggestions 
follow. 

What  affect  does  the  number  of  loos  yarded  p  e  r  ; y : ’ e 
have  on  productivity?  As  the  nymDer  z(  logs  yarded  ;er 
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I 

cycle  Increases,  hooking  time,  lateral  Inhaul  time,  and 
j  Inhaul  time  would  be  expected  to  Increase.  The  total 

number  of  cycles  needed  to  complete  the  operation  however, 
would  decrease.  How  many  logs  should  be  yarded  per  turn? 

I  What  effect  does  the  type  of  cable  logging  machine 

used  have  on  the  total  cost  of  the  operation?  For  example, 
a  very  powerful  machine  would  be  expected  to  shorten  Inhaul 
j  and  lateral  Inhaul  time.  Perhaps  more  logs  could  be  yarded 

per  turn.  However,  the  hourly  machine  cost  would  be 
expected  to  Increase  as  the  machine  capacity  Increases. 
Would  the  total  harvesting  time  for  the  entire  operation 
decrease  enough  to  warrant  the  use  of  this  larger  machine? 
In  addition  to  modeling  these  machine  specific 
^  variables,  a  subroutine  to  calculate  the  cost  of  the 

operation  Is  needed.  In  general,  the  least  expensive 
system  alternative  will  be  chosen. 

I 

Tho  size  of  the  labor  force  should  be  considered  when 
predicting  Independent  service  times.  What  is  the  effect 
of  having  a  larger  labor  force?  Does  the  total  cost  of  the 

I 

operation  decrease? 


Site  specific  variables  such  as  slope,  brush  1  ness,  and 
stand  structure  should  be  considered.  In  general,  a  site 
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A  subroutine  to  build  a  stand  of  trees  to  be  harvested 
Is  needed.  Cartesian  coordinates,  a  diameter,  height, 
and  weight  for  each  tree  could  be  estimated  from  cruise 
data.  A  decision  rule,  specified  by  the  user,  would  then  be 
employed  to  determine  which  tree  or  trees  should  be  yarded 
on  a  given  turn.  This  would  add  to  the  realistic  aspect  of 
the  operation  and  many  additional  "what-lf"  questions  could 
be  addressed. 

For  example,  many  different  logging  configurations 
could  be  considered.  Several  feasible  log  landing 
locations  could  be  examined  and  the  best  one  chosen.  The 
question  of  how  closely  the  corridors  should  be  spaced  as 
opposed  to  Increasing  lateral  outhaul  distances  could  be 
examined  In  greater  detail  than  addressed  In  Appendix  B. 

With  some  of  these  suggestions  for  further  research 
Implemented,  the  simulation  model  has  the  potential  to 
become  a  valuable  tool  for  planning  cable  logging 
operations  In  the  Allegheny  Region. 


LITERATURE  CITED 


Amerl can-Pul pwood  Association,  Harvesting  Research  Project. 
1972.  User  guide  for  the  Harvesting  System  Simulator, 
production  and  operations  report,  cost  analysis  report 
programs.  Atlanta,  Georgia. 

Bare,  B.  B.,  B.  A.  Jayne,  and  B.  F.  Anholt.  1976.  A 
simulation-based  approach  for  evaluating  logging 
residue  handling  systems.  USDA  Forest  Service  General 
Technical  Report  PNW-45. 

Baumgras,  J.  E.,  and  P.  A.  Peters.  1985.  Cost  and 

production  analysis  of  the  Bltteroot  mlnlycrder  on  an 
Appalachian  hardwood  site.  USDA  Forest  Service  Research 
Paper  NE-557. 

Biller,  C.  J.,  and  E.  L.  Fisher.  1984.  Whole-tree 
harvesting  with  medium  capacity  cable  yarder. 
Transactions  ASAE,  Vol  .  27(1),  pp.  2. 

•Biller,  C.  J.,  and  P.  A.  Peters.  1982.  Testing  a 

prototype  cable  yarder  for  tree  thinning.  Transactions 
ASAE,  Vol.  25(4),  pp.  901. 

Biller,  C.  J.,  and  P.  A.  Peters.  1984.  Harvesting  whole- 
tree  Appalachian  hardwoods  with  a  Washington  78  yarder. 
Transactions  ASAE,  Paper  No.  83-1605. 

Bradley,  D,  P.,  F.  E.  Blltonen,  and  S.A.  Wlnsauer.  1976. 

A  computer  simulation  on  full-tree  field  chipping  and 
trucking.  USDA  Forest  Service  Research  Paper  NC-129. 

Bradley,  D.  P.,  and  S.  A.  Wlnsauer.  1976.  Solving  wood- 
chip  transport  problems.  USDA  Forest  Service  Research 
Paper  NC-138. 

Bussel,  W.  H.,  J.  N.  Hool,  A.  M.  Leppert,  and  G.  R. 

Harmon.  1969.  Pulpwood  harvesting  systems  analysis. 
Interim  report  to  the  Southern  Executives  Association. 
Auburn  University,  Auburn,  Alabama. 

Cabbage,  F.  W.,  and  A.  H.  Gorse.  1975.  Mountain  logging 
with  the  Bltteroot  Mlnlyarder.  Mountain  Logging 
Symposium  Proceedings.  West  Virginia  University, 
Morgantown,  West  Virginia,  pp.  81. 

Fisher,  E.  1.,  H.  G.  Gibson,  and  C.  J.  Biller.  1980. 

Production  and  costs  of  a  live  skyline  yarder  tested  In 
Appalachia.  USDA  Forest  Service  Research  Paper  NE-465  . 


Fisher*  E.  L.*  D.  L.  Gouchenour,  and  C.  J.  Biller.  1984. 
Significant  factors  affecting  performance  of  a  Urus 
cable  yarder.  Transactions  ASAE*  Vol .  27(4),  pp.  962. 

Goulet*  D.  V.,  R.  H.  Iff,  and  D.  L.  Slrols.  1979.  A 

survey  of  timber  harvest  simulation  models  for  use  In 
the  south.  USDA  Forest  Service  General  Technical  Report 
SO-25  . 

Goulet,  D.  V.,  R.  H.  Iff,  and  D.  L.  Slrols.  1980. 

Analysis  of  five  forest  harvesting  simulation  models, 
part  II:  Paths,  pitfalls,  and  other  considerations. 
Forest  Prod.  J.  30(8):18-22. 

Hines,  G.  S.,  C.  Y.  Liu,  D .  B.  Webster,  and  D.  L.  Sirois. 
1983.  Features  of  a  skidding  module  In  a  timber 
harvesting  computer  simulation  model.  For.  Prod.  J. 
33(7/8) : 11-16. 

Hines,  G.  S.,  D.  B.  Webster,  D.  V.  Goulet,  and  D.  L. 

Slrols.  1981.  A  second  generation  of  timber 
harvesting.  Simulation.  Nov.  1981:165-168. 

Hocklnson,  J.  1986.  Timber  management  staff,  Allegheny 
National  Forest,  Warren,  Pennsylvania.  Personal 
communication.  11-13-86. 

Johnson,  L.  R.  1970.  Simulation  of  the  loading  and  hauling 
subsytems  of  a  logging  system.  M.  S.  Thesis, 

Dept,  of  Forestry,  Montana  State  Univ.,  Bozeman, 

Montana.  235  p. 

Johnson,  L.  R.,  C.  J.  Biller,  and  D.  L.  Gouchenour,  Jr. 

1972.  Simulation  analysis  of  t i mb e r- h a r vest  1 ng  systems. 
Proceedings  of  the  23rd  Annual  AIIE  Conference,  pp.  353. 

Johnson,  L.  R.,  and  C.  J.  Biller.  1973.  Wood-chipping 

and  a  balanced  logging  system:  Simulation  can  check  the 
combinations.  Winter  Meeting  ASAE,  Paper  No.  73-1537. 

Johnson,  L.  R.  1976.  SAPLOS  documentation  and  use. 

Final  Report  to  USDA  Forest  Service,  NE  Forest 
Experiment  Station.  Univ.  of  Idaho,  College  of 
Forestry,  Wildlife  and  Range  Science,  Moskow,  Idaho. 

Killham,  J.  R.  1975.  The  development  of  a  forest 

harvesting  simulation  model.  M.  S.  Thesis,  Dept,  of 
Industrial  Engineering,  Auburn  Univ.,  Auburn, 

Alabama.  158  p. 

Kleindorfer,  G.  B.  1986.  Management  Science  Dept,  of  The 
Pennsylvania  State  University,  University  Park,  PA. 
Personal  communication.  5-7-86. 


Law,  A.  M.,  and  W.  D.  Kelton.  1982.  Simulation  modeling 
and  analysis.  McGraw-Hill,  Inc.  New  York,  NY. 

LeDoux,  C.  B.  1985.  When  Is  hardwood  cable  logging 
economical?  Journal  of  Forestry,  May  1985. 

LeDoux,  C.  B.,  and  D.  A.  Butler.  1981.  Simulated  cable 
thinning  In  young-growth  stands.  Forest  Scl.  27(4): 
745-757. 

LeDoux,  C.  B.,  and  L.  W.  Starnes.  1986.  Cable  logging 

production  equations  for  thinning  young-growth  Douglas- 
f 1 r .  Forest  Prod.  J.  36 ( 5 ) : 2 1-24 . 

Liu,  C.  Y.  1981.  Development  of  a  skidding  module  for  a 
timber  harvesting  simulation  model.  Unpublished  M.  S. 
thesis.  Department  of  Industrial  Engineering,  Auburn 
Unlv.,  Auburn,  Alabama.  189  p. 

Lord,  R.  G.  1985.  A  Forplan  model  for  aggregate  timber 
supply  analysis  In  the  Allegheny  Region  of 
Pennsylvania.  Unpublished  M.  S.  Thesis,  School  of 
Forest  Resources,  The  Pennsylvania  State  Unlv., 
University  Park,  Pennsylvania.  265  p. 

Martin,  A.  J.  1975.  THATS:  Timber  and  Harvest  Transport 
Simulator.  USDA  Forest  Service  Research  Paper  NE-316. 

McGaughey,  R.  J.  1983.  Microcomputer  simulation  of  skyline 
logging  systems.  Unpublished  M.S.  thesis.  Department  of 
Forestry,  Purdue  Unlv.,  West  Lafayette,  Indiana. 

212  p. 

Naylor,  T.  H.,  and  J.  M.  Finger.  1967.  Verification  of 
computer  simulation  models.  Management  Science 
14(2) :B92-B101. 

O'Hearn,  S.  1977.  Economic  and  performance  comparison 

between  full  tree  chipping  and  conventional  harvesting 
systems  on  a  variety  of  stand  types  In  the  southeast. 
M.S.  Thesis.  Dept,  of  Forestry,  Va.  Polytech.  Inst, 
and  State  Unlv.,  Blacksburg,  Virginia.  128  p. 

Padgett,  M.  L.  1982.  A  statistical  methodology  for  the 
development  of  computer  simulation  models  for  forest 
harvesting.  Unpublished  M.  S.  Thesis,  Dept,  of 
Industrial  Engineering,  Auburn  University,  Auburn,  AL. 
546  p  . 

Pegden,  C.  D.  1985.  Introduction  to  SIMAN  with  version  3.0 
enhancements.  Systems  Modeling  Corporation,  State 
College,  Pennsylvania. 


55 


j 


•a 


i 


5 


\ 

r 


Peters,  P.  A.  1984.  Steep  slope  clearcut  harvesting  with 
cable  yarders.  Proceedings  of  Harvesting  Small  Trees  In 
the  South,  April  18-20,  BIT  1x1,  Mississippi. 

Peters,  P.  A.  1985.  Cost-estimating  programs  for  cable 
yarding.  Winter  Meeting  ASAE,  Paper  No.  85-1590. 

Peters,  P.  A.,  and  J.  E.  Baumgras.  1984.  Production  and 
cost  analysis  of  the  Bltteroot  mlnlyarder  In 
Appalachia.  Winter  Meeting  ASAE,  Paper  No.  84-1603. 

Powell,  D.  S.,  and  T.  J.  Consldlne.  1982.  An  analysis  of 
Pennsylvania's  forest  resources.  USDA  Forest  Service 
Resource  Bulletin  NE-89. 

Prltsker,  A.  A.  B.  1974.  The  GASP  IV  simulation  language. 
John  Wiley  and  Sons,  Inc.,  New  York. 

Prltsker,  A.  A.  B.  1984.  Introduction  to  simulation  and 
SLAM  II.  John  Wiley  and  Sons,  Inc.,  New  York. 

Prltsker,  A.  A.  B.,  and  P.  J.  Kavlat.  1964.  Simulation 
with  GASP  II:  a  fortran  based  simulation  language. 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  New  Jersey. 

Rodgers,  J.  H.  1984.  An  Input  data  p re- p rocesso r  for  a 

timber  harvesting  computer  simulation.  Design  Project, 
Dept,  of  Industrial  Engineering,  Auburn  Unlv.,  Auburn, 
AL.  163  p . 

Rossle.  K.  M.  1983.  A  case  study  of  the  Koller  K-300 
yarder  on  a  national  forest  timber  sale  In  the 
Appalachian  region.  M.S.  Thesis.  Dept,  of 
Forestry,  Va.  Polytech.  Inst,  and  State  Unlv., 
Blacksburg,  Virginia.  115  p. 

Rummer,  R.  B.  1986.  Graduate  Co-op  Research  Engineer. 

USDA  Forest  Service,  Southeast  Forest  Experiment 
Station,  Auburn,  Alabama.  Personal  comun Icat 1  on . 
7-10-86 . 

Sargent,  R.  G.  1984.  A  tutorial  on  verification  and 

validation  of  simulation  models.  Proceedings  of  the 
IEEE  Winter  Simulation  Conference,  pp.  115. 

Schrlber,  T.  1974.  Simulation  using  GPSS.  John  Wiley  and 
Sons,  Inc.,  New  York. 

Shannon,  R.  E.  1975.  Systems  simulation:  the  art  and 

science.  Prentice-Hall,  Inc.,  Englewood  Cliffs,  New 
Jersey. 

Shruben ,  l.  W,  1980.  Establishing  the  credibility  of 
simulations.  Simulation  38  ( 2  )  :  10  1- 1  05  . 


•.Am.  w.*— 


.  -  V  /  *  \ 


.  -  ■  •  .  v  v  „■% 
■-u  *  -  ”  — 11  *  A  -  '■ 


•fcwwLih 


>  a  *  mwAwJ  m  »'i  i  *,'■  w. 


Stark,  J.  I.  1968.  A  Simulation  Model  for  the  Common 
Pulpwood  Harvesting  Systems  of  the  Southern  Pine 
Region.  M.  S.  thesis.  Dept,  of  Industrial 
Engineering,  Georgia  Inst,  of  Tech.,  Atlanta,  Goergla. 
196  p. 

STSC,  Inc.  1985.  Statgraphlcs  statistical  graphics 
system.  Rockville,  MD . 

Webster,  D.  B.  1975.  Development  of  a  flexible  timber 

harvest  simulation  model.  Forest  Prod.  J.  25(15:40-45. 

Webster,  D.  B.»  M.  L.  Padgett,  G.  S.  Hines,  and  D.  L. 

Slrols.  1984.  Determining  the  level  of  detail  In  a 
simulation  model-a  case  study.  Comput.  <5,  Ind.  Engng. 
8(3/45:215-225. 

Webster,  D.  B.,  M.  L.  Padgett,  and  D.  L.  Slrols.  1983. 

Features  of  a  felling  module  using  a  f el  1 e r-b u nche r  In 
a  timber  harvesting  computer  simulation  model.  For. 
Prod.  J.  33(65:11-16. 


The  Forest  Harvesting  Simulation  Model  (FHSM)  is  a 
FORTRAN  /  GASP  II  (Prltsker  and  Kavlat  1964),  event- 
oriented  model  that  can  be  used  to  simulate  a  wide  variety 
of  saw  log  and  pulpwood  harvesting  operations  In  the  South 
(Klllham  1975).  Webster  (1975)  reported  that  the 
requirements  In  developing  this  model  were  that  It  (a)  be 
flexible  enough  to  duplicate  the  major  systems  used  In 
timber  harvesting  In  the  South,  (b)  be  detailed  enough 
allow  for  analysis  of  Individual  harvesting  functions,  and 
(3)  possess  a  high  "degree  of  bel 1 evab 11 1 ty"  In  the  way  It 
duplicates  a  system’s  operations. 

Webster  further  describes  the  model.  To  satisfy  the 
first  objective,  the  model  simulates  the  functional 
elements  such  as  felling,  limbing,  bucking,  skidding,  etc. 
In  various  configurations.  While  the  configurations  do  not 
cover  all  types  of  harvesting  systems  In  the  South,  they  do 
cover  a  wide  range  of  them. 

Each  function  has  been  built  as  a  'separate  component 
In  the  model.  Development  of  a  model  formulation  for  a 
specific  system  Involves  assembling  the  various  components 
that  comprise  that  system. 

To  satisfy  the  second  objective,  components  of  the 
model  (the  harvesting  functions)  are  defined  so  as  to  allow 
different  pieces  of  equipment  to  perform  the  same  function, 
but  also  differentiates  the  equipment  characteristics  and 


work  capabilities.  For  example,  two  different  types  of 
skldders  may  be  used  In  the  skidding  function. 


To  satisfy  the  third  objective,  the  model  Is  detailed 
enough  to  allow  for  the  user  to  follow  the  flow  of  wood 
through  the  model.  In  order  to  achieve  this,  a  tree  Is 
felled  and  It  provides  the  basis  for  wood  Input  to  the  rest 
of  the  model . 

Model  Input  can  be  divided  up  Into  general  and 
operational  categories.  General  Input  Influences  the 
performance  of  all  of  the  harvesting  functions  (e.g.  tree 
stand  mix,  tree  size,  merchantable  height,  etc.). 
Operational  Input  Is  that  which  Influences  one  harvesting 
function  (e.g.  skldder  capacity,  travel  rate  of  skidder, 
etc.  ) . 

In  general,  the  model  output  consists  of  production 
data  In  board  feet,  cubic  feet,  and  weight  for  each 
harvesting  component  over  the  time  horizon  simulated. 
Additionally,  for  each  harvesting  function,  production 
statistics,  productive  time,  and  Idle  time  are  kept  on  each 
piece  of  equipment  or  crew  member. 

Ful  1  -T ree  Field  Chipping  and.  T ransport  $ jJiuljjLar 

The  Full-Tree  Field  Chipping  and  Transport  Simulator 
(FTFC)  Is  composed  of  two  GPSS  (Schrlber  1974)  simulation 
models,  one  for  chipping  (Bradley  et  al.  1976)  and  one  for 
transport  (Bradley  and  Wlnsauer  1976),  and  was  designed  to 
duplicate  the  features  of  a  full-tree  field  chipping 
operation  Including  a  stand  of  trees,  f e 1 1 e r- b u nc he r s , 


skldders,  a  chipper  with  loader,  trucks  and  vans,  one  or 
two  optional  setout  trucks,  a  mlllyard  scale  and  a  chip 
dumper. 

Bradley  et  al .  (1976)  provide  a  model  description. 

The  model  attempts  to  mirror  the  harvesting  operation  down 
to  elemental  time  within  each  function  In  order  to  give  the 
model  a  high  degree  of  Integrity.  It  Is  machine 
Independent  In  that  machine  speeds  and  capacities  may  be 
altered  to  account  for  different  equipment  brands  and 
sizes. 

FTFC  does  not  consider  delays  caused  by  machine 
breakdown  but  has  Instead  concentrated  on  delays  caused  by 
machine  Interaction.  For  example,  the  skldder  may  "look 
ahead"  to  see  If  another  skldder  Is  unloading  at  the 
chipper.  If  not.  It  will  take  Its  load  to  the  chipper.  If 
another  skldder  Is  unloading,  however.  It  will  deposit  Its 
load  at  a  stockpile  and  bring  It  to  the  chipper  later,  thus 
representing  a  system  delay. 

Wood  flows  through  FTFC  In  a  manner  Identical  to  FHSM 
with  one  exception.  Instead  of  actually  passing  wood  from 
one  function  to  the  other,  parallel  random  number  streams 
are  generated  at  each  function.  Wood  flow  Is  Identical  to 
that  of  FHSM  when  wood  passes  through  each  function  In  the 
same  sequence.  However,  when  the  orders  of  two  bunches 
have  been  Inverted  from  felling  to  skidding,  the  trees  In 
the  bunches  have  been  effectively  rearranged,  and  hence, 
skid  turn  statistics  will  vary  from  model  to  model. 


Model  Input  requirements  Include  tree  locations  (In  x, 
coordinates)*  tree  volumes,  and  felling  order  of  each  tree 
In  the  stand.  Elemental  machine  productivity  and  capacity 
data  are  also  required. 

The  model's  report  generator  provides  detailed  production 
and  cost  statistics  by  operation,  system  energy 
consumption,  and  net  energy  produced,  the  latter  two  In  the 
form  of  BTU's  (Goulet  et  al  .  1979  ). 

Harvest  System  Simulator 

The  Harvest  System  Simulator  (HSS)  Is  a  FORTRAN-based 
time  and  event  oriented  simulation  program  and  Is  part  of  a 
larger  package  known  as  the  Harvesting  Analysis  Technique 
(HAT)  (American  Pulpwood  Association  1972).  The  model  was 
designed  to  simulate  the  productive  and  non-productive 
(dow  Ime,  breaks,  etc.)  activities  of  a  harvesting  system 
by  simulating  the  Interaction  between  harvesting  equipment 
and  the  stand  being  harvested. 

HSS  Is  different  from  FTFC  and  FHSM  in  that  it  focuses 
on  the  larger  systems  aspect  of  the  harvesting  operation, 
l.e.,  function  to  function  Interaction  and  not  the  detail 
of  any  one  function  (Goulet  et  al.  1980). 

O'Hearn  (1977)  described  HSS.  A  maximum  of  14 
machines,  working  In  any  combination,  and  a  maximum  of  six 
aggregations  of  like  machines  (phase)  can  be  simulated. 

The  harvested  tract  can  be  divided  into  a  maximum  of  14 


harvesting  areas  that  can  differ  In  stand  type,  volume  per 
acre,  species,  composition,  and  skidding  distance  to  the 


primary  landing.  Individual  harvesting  areas  have  no 
acreage  or  volume  limits.  Unique  production  rates  may  be 
specified  for  each  harvesting  area-machine  combination. 

The  user  controls  the  order  In  which  harvesting  areas  are 
processed.  Terrain  and  stand  limitations  are  modeled 
through  move  or  travel  rate  modifiers  and  deck  locations. 
Wood  flows  from  phase  to  phase  In  aggregated  volumes. 

Nonproductive  activities,  such  as  machine  failures, 
breakdowns,  delays,  and  servicing,  can  be  Imposed  logically 
or  stochastically.  Repairs  can  be  made  at  the  stump,  the 
deck,  or  the  shop.  Repairs  at  the  stump  hold  the  machine 
In  place,  while  repairs  at  the  deck  or  shop  require  the 
machine  to  move  to  the  primary  deck.  Delays  are  divided 
Into  two  types:  major  and  minor.  Major  delays  bring  the 
machine  back  to  the  deck,  while  minor  delays  leave  It  In 
place.  A  distribution  between  productive  and  nonproductive 
time  can  be  provided. 

Program  output  Is  provided  by  two  report  generators 
that  can  be  called  on  separately  or  jointly.  They  provide 
time,  production,  cost,  and  revenue  statistics.  Also, 
discounted  cash  flow  and  return  on  Investment  analysis  can 
be  made.  All  output  reports  and  detailed  and  complete. 

Slmul  at  1  o n  App.l.l _e4  la  Lamina  5. y -terns 
Simulation  Applied  to  Logging  Systems  (SAPLOS) 
(Johnson  1976)  Is  a  discrete  event,  FORTRAN  /  GASP  IV 


(Prltsker  1974),  general  harvest  simulation  model  that  Is 
adaptable  to  a  number  of  harvesting  systems. 


Johnson  et  a  1 .  (1972)  reported,  the  development  o ^  tne 

1 1  mb e r- h a r v e s 1 1 n g  model  orogressed  through  three  phases  of 
analysis:  (1)  Identify  and  classify  the  most  common 

logging  systems  in  Appalachia  In  terms  of  the  general 
operations  (subsystems  themselves)  Involved  and  their 
points  of  Interaction.  (2)  Within  each  subsystem  Identify 
the  activities  particular  to  that  subsystem  adn  determine 
the  sequence  of  these  activities.  Activities  are  the 
various  operations  of  a  subsystem  such  as  movement  of  a 
skldder  to  the  landing  or  the  hook-up  of  logs  by  the 
skldder.  (3)  Combine  the  particular  activities  of  the 
subsystems  Into  the  general  events  and  activities  Involved 
In  the  actual  computer  model.  Events  signal  the  beginning 
and  ending  of  activities. 

Johnson  further  describes  the  model.  In  the  first 
phase  of  development,  six  subsystems  were  Identified  and 
represents  the  standard  operations  of  this  model.  They  are 
felling,  bucking,  prebunching,  skidding,  loading,  and 
hauling. 

In  the  model,  the  operation  of  a  particular  subsystem 
Is  represented  by  the  activities  that  make  up  that  operation 
and  the  events  that  signal  activity  changes.  The  second 
phase  of  the  model  development  was  accomplished  through 
the  1 ndent 1 f 1  cat  1  on  of  these  activities.  For  example,  the 
arrival  of  a  transporting  vehicle  at  either  the  stump  or 
skldroad  signals  the  first  skidding  activity-winching  or 
hooking  the  log  to  the  sk 1 dder.  The  stump  depicts  the  area 
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wnere  fie  tree  was  cut.  An  e n c - o f - h o o *  1  n g  event  marks  the 
end  of  the  « inching  activity  and  the  beginning  of  a  move  to 
either  another  tree  or  the  landing,  depending  on  whether  or 
not  the  load  Is  full. 

The  third  phase  of  development  was  accomplished 
throuqn  the  Identification  of  five  "cr'tlcal  locations" 
where  the  si*  logging  operations  overlap.  A  "critical 
location"  on  the  ground  forms  a  control  point  in  the 
simulation  model.  They  are  (1)  the  stump  or  tree  location 
In  the  woods,  (2)  the  skidroad  (when  precunchers  are  used), 


( 3 : 

the  landing. 

(  4  ) 

tne 

prenaul  cock 

(  when 

shuttle  truck 

s 

and 

prehaul  trail 

e  r  s 

are 

used  and,  (5! 

the 

process  1 ng  po  * 

n  t 

such  as  a  mill  yard). 

Two  stages  of  input  are  required  fn  the  computer 
model.  The  two  stages  consist  of  objective  information  and 
subjective  information.  Objective  Information  consists  of 
the  distributions  used  in  the  production  equations  and  are 
used  to  obtain  travel  distances  and  load  sizes.  Subjective 
Information  Is  used  In  equations  to  calculate  production 
times.  The  subjective  Information  describes  the  ’egging 
system  being  simulated  In  terms  of  the  variables  un;que  to 
that  system,  e .  g .  ,  a  description  of  the  terra 'in  (  n  terms  of 
slope  and  soil. 

Model  output  consists  of  jroduct  ’  on  f  r  q  m  the  operation 
In  cubic  feet  of  timber,  t  n  e  time  required  to  achieve  t h i s 
production,  and  the  total  cost  involved  4  *  pro  due  •  o  q  the 


timber, 


t imber  j±arv£ii  ana  iraiianaxi  Him J .alar 

The  Timber  Harvest  and  Transport  Simulator  (THATS) 
(Martin  1975)  Is  a  FORTRAN  based  time  and  event  oriented 
simulation  program  which  models  the  major  harvesting 
systems  of  the  Appalachian  Region.  Martin  (1975)  describes 
the  structures,  methodology,  and  main  components  of  THATS. 
THATS  Is  built  around  a  main  program  composed  of  eight 
components  (felling,  bunching,  skidding,  bucking,  loading, 
hauling,  roadbulldlng,  and  cost  accounting)  and  a  "clock." 
The  model  Is  a  time  oriented  simulation  In  which  simulated 
time  on  the  clock  Is  advanced  one  minute,  then  checks  are 
made  for  active  events. 

Simulated  event  times  are  generated  either  from  gSen 
averages  and  standard  deviations,  or  from  event  times 
oroduced  by  a  regression  equation  developed  from  collected 
data.  All  random  variable  event  times  have  either  a  norma’ 
or  a  lognormal  distribution.  If  any  shewing  is  present  in 
tie  t'me  study  data,  the  lognormal  distribution  j  s  used. 
Regression  equations  from  data  collected  for  Appal  ach'  an 
logging  operat  ''ons  a  r  e  contained  1  n  the  report. 

The  system  simulates  one  day  at  a  t  * m e  am!  s  *  u  t  s  :  :  w n 
at  the  end  of  t  me  w  o  r i  n  g  day  in  a  stay  cere  d  manner,  w;tK 
e  a  r  r  e w  f  '  n  i  s  n  i  n  ;  tne  day  c  ’  c  s  e  to  g  u  ‘  1 1  •  n  ;  ♦  '  me  ,  t  n  o  j  ;  * 
d  e o  e r  y  "  c  on  the  t  a  s  ■  at  *  a  n  d  ,  seme  ~ a  ,  r  °  a  ■'  1 1  ’  e  ear', 
ant  some  m  a  /  be  a  ’  1 1 ' e  'ate. 


and  the  output  from  the  operation  Is  a  volume.  The  next 
operation  draws  trees  or  pieces  from  the  deposited  output 
volume,  but  these  new  trees  or  pieces  have  no  relation  to 
the  Input  trees  of  the  first  operation  except  that  their 
total  volume  equals  the  doposlted  output  volume.  For 
example,  a  tree  Is  generated  for  the  felling  operation, 
felling  statistics  are  collected  on  that  tree,  and  the 
volume  of  that  tree  Is  deposited  for  the  bunching 
operation.  The  bunching  operation  will  now  generate  new 
trees  up  to  the  volume  of  the  felled  trees  deposited. 

Input  for  THATS  Includes  time  and  motion  data 
appropriate  to  teh  system  being  studied.  For  those 
functions  In  which  operating  or  delay  time  has  not  been 
related  quantitatively  to  system  or  stand  parameters,  the 
Input  consist  of  expected  values  (and  their  standard 
deviations) . 

Output  from  THATS  Include  system  status  Information, 
time  summaries,  production  summaries,  and  cost  summaries 


Martin  ( 1975 ) . 


Comoar 1  son  q±  Two  Harvest  Unit  Configurations 

The  simulation  model  was  used  to  compare  two  harvest 
unit  con f 1 gu rat  1  on s .  The  utility  of  the  model  Is 
highlighted  and  the  procedure  for  modifying  the 
experimental  frame  Is  Illustrated. 

In  the  first  scenario,  the  unit  was  harvested  from  four 
corridors  at  landing  1  and  two  corridors  at  landing  2;  in 
the  second  scenario  there  were  two  corridors  at  each 
landing.  Minimization  of  make-span  (total  time  to  complete 
the  operation)  was  used  as  the  performance  criteria. 

The  crux  of  the  analysles  was  determining  whether  <t 
takes  more  time  to  harvest  the  un ‘t  from  more  corridors  Dut 
with  relatively  short  lateral  hauls  versus  harvesting  from 
fewer  corridors  but  with  relatively  long  lateral  hauls. 

The  first  scenario  was  modeled  us'nq  equations 
developed  from  the  time  and  motion  study  that  was 
carried  out.  Modifications  were  then  made  to  three  c* 
these  equations  In  order  to  reflect  the  changes  that 
made  for  the  second  scenario. 

First,  the  times  (number  of  cycles  to  complete  at  a 
given  o  u  t  h.  a  u  1  distance)  empirical  distr'oution  was 
modified  to  reflect  the  increased  1  i  x  1  1  h  o  o  d  of  hav<r:;  me.  *-  e 
yarding  cycles  at  a  q'ven  distance  than  ' n  the  first 
scenar'o.  The  lateral  Inhaul  and  1 atera’  c  u  t  r  a  u  '  equal 
were  then  modified  to  reflect  the  increased  '  1  > '  ■ hoi  •  o * 
having  relatively  longer  lateral  hau's.  T  n  P  u  P  ■-  K  a r  ;  e  •  a  ’ 


illustrated  Dy  compar’ng  the  two  r,  ( v  a  g  p.»ppr<rrPnta 


r  ^  fr- 


listings.  The  listing  for  scenario  1  (pp.  73-74)  Is 
presented  In  Its  entirety.  The  listing  for  scenario  2 
(p.  75)  only  Includes  the  changes  that  were  made.  All 
other  code  Is  Identical  to  that  of  the  first  listing.  The 
SIMAN  model  frame  listing  (Identical  for  both  scenarios) 

Is  presented  (pp.  76-79).  Directions  for  the  modification 
of  the  experimental  frame  are  presented  below. 

Interpretation  Model  Qylflyj 
When  Input  to  a  simulation  model  are  random  (as  Is  the 
case  here)  the  output  statistics  will  vary  from  run  to  run. 
It  Is  thus  necessary  to  perform  several  replications  of 
each  scenario  and  perform  statistical  analyses  on  the  data 
of  'nterest.  The  procedure  Is  straightforward  and  will  not 
he  presented  here.  The  output  from  two  simulation  runs 
(one  from  each  scenario)  are  presented  In  Table  B.l  and  B.2 
and  are  sufficient  to  highlight  the  utility  of  the  model. 

* r om  these  tables  It  can  be  seen  that  make-span  for 
scenar'o  one  was  144,400  seconds  (40.1  hours)  and  for 
scenario  two  127,000  seconds  (35.3  hours).  Therefore, 
scenar'o  two  would  be  chosen. 

Lire*  am  Far  M£.fl_L.LUAlian  qj_  LitsrlJUfiDiAl  Exams 
The  current  model  Is  set  up  to  model  a  harvesting 
colouration  w'th  two  landings.  Modification  of  the  SIMAN 
mode1  frame  would  be  necessary  to  consider  additional 
1  a  n  d  '  r  g  s  .  ^  n  i  s  d  '  s  c  u  s  s  1  o  n  w  i  1  1  thus  only  prov'de 

'  n  s  t  •"  j  r.  t  '  o  n  s  ‘or  the  modification  of  the  number  of 
o  r  r  '  ■’  o  r  s  -ar  vested  and  of  parameters  associated  w'th  those 


Table  B.l:  SI HAN  summary  statistics  for  scenario 
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corridors  (corridor  length,  starting  distance,  and 

Incremental  corridor  Increase).  In  general,  n  corridors 

may  be  considered  as  follows: 

1)  In  the  COUNTERS  element  change  the  last  number  to  n+l; 

2)  In  the  PARAMETERS  element  2,  change  the  second  number 
to  n  ; 

3)  In  PARAMETERS  element  5,  specify  the  Incremental 
corridor  Increase  desired; 

4)  In  PARAMETERS  element  6  and  11,  specify  the  last 
corridor  number  at  landing  1  and  n  respectively; 

5)  In  PARAMETERS  element  27  through  27+n,  specify  the 
corridor  lengths; 

6)  In  PARAMETERS  element  28+n  through  28+2n,  specify  the 
starting  distance  for  harvesting  at  each  corridor. 


##SIMAN  Experimental  Frame  -  Scenario  1** 

BEGIN; 

PRMECT,CAELE  YARCER,  PETER  HCCES, 2/ 12/1987; 

DISCRETE,5,3; 

TALLIES: 

1,  OUTHALL  DISTANCE,  11: 

2,  OUTWVJL  TINE,12: 

3,  LAT.  CUTH_.  TINE,  13s 

4,  KXK  TINE,  14: 


5,  LAT.  INR.  TDE,15: 

6,  IfmjL  TINE, 16: 

7,  INH_.  DLY.  TINE,  17: 

8,  UNKXKINS  TBE,18: 

9,  CLR.  LAND  TINE,  19; 
COJNTERS:l,END  CF  RJN,7; 
REPLICATE,  1; 

PARANETERS: 


1,2.0: 

2,6.0: 

3,14400.0: 

4,5400.0: 

5,25: 

6,4.0: 

7,4.46: 

8,0.05204: 

9,5.14: 

10,0.26617: 

11,6.0: 

12,  .925: 

13,100.63: 

14,56.59: 

15,2.96: 

16,634: 

17,0.03038: 

18,24.68: 

19,0.10185: 


I  TOTAL  NLNEER  CF  LANDINGS 
1  TOTAL  NLNEER  CF  CCERIDCRS 
I  LANDING  SETLP  TINE  (SBCCNDS) 

1  CCRRIDCR  SETIP  TINE  (SBCCNDS) 

I  INCREMENTAL  CXRRIDCR  INCREASE 
I  LAST  0CFS3DCR  MNEER  AT  LANDING  MNEER  1 
I  EETAO  FCR  LANDING  NLNEER  1 
I  EETA1  FCR  LANDING  NLNEER  1 
I  EETOO  FCR  LANDING  NLNEER  1 
I  BETH  FCR  LANDING  NLNEER  1 
1  LAST  CCRRIDCR  NLNEER  AT  LANDING  NLNEER  2 
I  PRCB.  CF  NO  INHALL  DELAY  AT  LANDING  NLNEER 
I  DURATION  CF  INWALL  DELAY  AT  LANDING  NLNEER 
!  LATERAL  CUT1NE  BETA  PARAM  (1)  FCR  LANDING 
!  LATERAL  CLHINE  ALFHA  PARAM  (2)  FCR  LANDING 
I  EETAO  FCR  LANDING  NLNEER  2 
I  EETA1  FCR  LANDING  NLNEER  2 
I  EETOO  FCR  LANDING  NLNEER  2 
I  EET11  FCR  LANDING  NLNEER  2 


20, .24,1, 34,2, .48,3, .55,4, .58,6, .69,7, .72,8, .76,9, .79,10, .82,11, 


.86,12,1,20:  1TTNES  DISTRIBUTION 

21,3530,2.56:  I  LATERAL  CUTWJL  TINE  DISTRIBUTION  INFCFWTICN 

22,53.81,1.9015:  !  HOCKING  TINE  DISTRIBUTION  INFORMATION 


23, .01,4, .1,11, .2,14, .3,17, .4,22, .5,30, .6,40, .7,68, 


.8,92, .9, 136,1. ,387: 

24..  892.0.0.1.0.94. 16: 
25,18.96,7.73: 

26. . 858.0.0.1.0.45.0: 
27,300.0: 

28,300.00: 

29,300.00: 

30,300.00: 


I  LATERAL  INWALL  TINE  DISTRIBUTION  INFO. 

I  INWALL  DELAY  TINE  INFORMATION 
!  UNLOCKING  TINE  DISTRIBUTION  INFCWATICN 
!  CLEAR  CF  LANDING  TINE  INFORMATION 
I  LENGTH  CF  CCRRIDCR  NLNEER  1 
!  LENGTH  CF  CCRRIDCR  MJNEER  2 
!  LENGTH  CF  CCRRIDCR  NLNEER  3 
!  LENGTH  CF  CCRRIDCR  NLNEER  4 
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31,300.00: 
32,300.00: 
20.00: 
20.00: 
20.00: 
20.00: 
20.00: 
20.00; 


I  LBGTH  CF  CCFRIDCR  hiJ^R  5 
1  UEHaTH  CF  CCFRIDCR  U*ER  6 
I  STARTDISTANCE  AT  CCFRIDCR  MJ€ER  1 
!  STARTDISTANCE  AT  CCFRIDCR  MJNEER  2 
I  STARTDISTANCE  AT  CCFRIDCR  NLMSR  3 
!  STAR1DISWCE  AT  CCFRIDCR  MM5ER  4 
!  STARTDI  STANCE  AT  CCFRIDCR  NLMER  5 
1  STARTDISTANCE  AT  CCRUDCR  6 


BEGIN; 


**SIMW  Experimental  Frame  -  Scenario  2** 


CONTERS:l,END  CF  RJN,5; 

PARAMETERS: 

2,*4.0:  1  TOTAL  MNEER  CF  CCFRIDCRS 

• 

6^2. 0:  1  LAST  CCFRIDCR  MNEER  AT  LANDING  NLNEER  I 

11,4.0:  1  LAST  CCRRIDCR  NJEER  AT  LANDING  MMBER  2 

• 

20, .21,1, .28,2, .39,3, .42,4, .48,6, .62,7, .68,8, .74,9, .77,80, .80,11, 
.86,12,1,20:  JUNES  DISTRIBUTION 

21,50.00,2.56:  1  LATERAL  OUTWJL  TINE  DISTRIBUTION  INFORMATION 


23, .01,4, .07,11, .14,14, .21,17, .35,22,  .49,30, .63,40, .77,68, 

.92,92, .97, 136,1., 387:  I  LATERAL  INmJL  TINE  DISTRIBUTION  INFO. 


27,300.0: 

28,300.00 

29,300.00 

30,300.00 

31,20.00: 

32,20.00: 

33,20.00: 

34,20.00; 


!  LENGTH  CF  CCRUDCR  NLMBER  1 
!  LENGTH  CF  CCFRIDCR  NLNEER  2 
I  LENGTH  CF  CCFRIDCR  MMBER  3 
I  LENGTH  CF  CCRRIDCR  NUMBER  4 
!  STARTDISTANCE  AT  OCFRIDCR  NJNEER  1 
!  STARTDISTANCE  AT  CCFRIDCR  NUMBER  2 
!  STARTDISTANCE  AT  OCRRIDCR  NM3ER  3 
STARTDISTANCE  AT  CCFRIDCR  NLMER  4 


***SIMW  Model  Fratre*** 


BEGIN; 

SYNONVMS: 

TDES  =  X(  1) : 
CEDE  =  X(2): 
INCRE  =  X(3): 
EETAO  =  X(4): 
EETA1  =  X(5) : 
EETOO  =  X(6) : 
EET11  =  X(7): 
LNU-B  =  X(8) : 
CNLNB  =  NC(1): 
MML  =  X(9): 
NLNBC  =  X(10): 
LTINE  =  X( 11) : 
CLENS  =  X(  12) : 
SDIST  =  X(  13 ) : 
TTINE  =  A(l): 
CDIST  =  X  ( 14 ) : 
TCLEA  =  X( 15 ) : 
OTINE  =  X(  16 ) : 
LOTTM  =  X(17): 
LITtlM  =  X(18): 
ITINE  =  X(19): 
HCJTIM  =  X(20): 
UHTTM  =  X(21) : 
HTTM  =  X(22): 
LCCRR  =  X(23) : 
LSUTI  =  X(24) : 
CSUH  =  X(25); 


!  NLNEER  CF  CTOLES  AT  A  SPECIFIED  OJimjL  DISTANCE 
!  NUMBER  CF  CVOES  CCM=LETED  AT  CURRENT  OJTVWJL  DISTANCE 
1  INCRENBNTAL  CCRRIDCR  DISTANCE  INCREASE 
I  Y-INTEFCEPT  IN  ECUATTCN  TO  PREDICT  CUTmJL  TINE 
!  CCEFF.  CF  CDIST  IN  ECL1ATICN  TO  PREDICT  OUTHAUL  TINE 
1  Y-INTEFCEFT  IN  ECLATICN  TO  FREDICT  INHAUL  TINE 
I  CCEFF.  CF  CDIST  IN  EOJATICN  TO  FREDICT  IN-WUL  TINE 
!  CURRENT  LANDING  NUMBER 
!  CURRENT  CCRRIDCR  NJNEER 
I  TOTAL  NLNEER  CF  UUDINGS 
!  TOTAL  NLNEER  CF  CCRRIDCRS 
!  ELAPSED  TINE  SINCE  CCRRIDCR  EEGAN 
I  LENGTH  CF  CURRENT  CCRRIDCR 

!  DISTANCE  WFEFE  HARVESTING  EEGAN  AT  CURRENT  CCRRIDCR 
!  MARKS  SMLATED  TINE  AT  KGINMING  CF  A  TURN 
!  CURRENT  QLfTHALL  DISTANCE 
I  TINE  TO  CLEAR  CURRENT  LANDING  CF  LOGS 
!  CLTTHALL  TINE 
1  LATERAL  OJTmjL  TINE 
l  LATERAL  INHAUL  TINE 
1  INHALE  TINE 
!  HOCKING  TINE 
I  UNHOCKING  TINE 
!  INHAUL  DELAY  TINE 

I  NUNEER  CF  LAST  CCRRIDCR  AT  TVE  CURRENT  LANDING 
!  LANDING  SETUP  TINE 
CCRRIDCR  SETUP  TINE 


START  CREATE,  1; 

ASSIGN:  'NUEL'  =  00(1) 
ASSIGN:  'NLM3C'  =  C0(2) 
ASSIGN:  'LSUH'  =  COG) 
ASSIGN:  'CSUH'  =  00(4) 
ASSIGN:  'INCRE'  =  COG) 


ASSIGN  NINEL 
ASSIGN  MJ-EC 
ASSIGN  LSUTI 

assign  csun 

ASSIGN  INCRE 


-  ■  »  *  j  j  V* 
IfaiAiA'jj  L.V  Jr  . 


XXXXMMXXXXXKXXXXX***X*XXXXXXXXXXXXXXXXXXXXXXX»tt*-X*XKXXXXXXXXXXXXKXXXXXXXX 

*  INITIALIZE  VALUES  FCR  LANDING  fUMEER  CNE  * 

iHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHt  *****  *********************** 


LANDI  ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 
ASSIGN 


'LNLNB'  =  1; 
'LCCRR'  =  00(6) 
'EETAO'  =  00(7) 
'EETA1'  =  00(8) 
'EETO0'  =  00(9) 


ASSIGN  LNLNE 
ASSIGN  LCCRR 
ASSIGN  EETAO 
ASSIGN  EETA1 
ASSIGN  EETOO 


•XX*  A‘‘‘W,k  M  mV.  K -V  -  »  -  -  uV  JV 


'-V-L^Vv'’' V-  -  V  .  •  .-U'.  V 

!i,ki>  iA.*?*  n.,ft  MiJfc  aJfckJL wTkifcX  m MuJjk** m.j, 


ASSIGN:  'BETH'  =  00(10);  ASSIGN  EET11 

DELAY:  'LSUTI':  LEXT  ( NEWCCRR ) ;  DELAY  BY  LCUTI  TEEN  PROCEDE  TO 

BLOCK  LEWCCRR 


LAND2 


ASSIGN:  'LCCRR'  =  CO(Il) 
ASSIGN:  P(24,l)  =  00(12) 
ASSIGN:  P(24»4)  =  00(13) 


ASSIGN:  P(21,l)  =00(14); 

ASSIGN:  P(21»2)  =00(15); 

ASSIGN:  'EETAO'  =00(16); 

ASSIGN:  'BETA1'  =  00(17); 

ASSIGN:  'BETOO'  =  00(18); 

ASSIGN:  'EET11'  =00(19); 

DELAY:  'LSUTI':  LEXT(CCNTINUl); 


ASSIGN  LCCRR 

ASSIGN  FR(B.  CF  AN  INHALL  DELAY 
ASSIGN  CURATTON  CF  AN  INHALL 
DELAY 

ASSIGN  LATERAL  CUTmjL  DISTRIBUTION: 

1.  EETA  PARAMETER 

2.  ALPHA  PARAMETER 
ASSIGN  EET/O 
ASSIGN  EETA1 
ASSIGN  EETOO 
ASSIGN  EET11 

DELAY  BY  LSUTI  TTEN  PROCEDE  TO 
BLOCK  C0NTIMJ1 


NEWCCRR  COUNT:  1,1;  INCREMENTS  CNLNB 

BRANCH,  1: 

IF,  'CMNB'  .Ed  ’LCCRR’  +  1,  LEWLAND: 

ELSE,  C0NTINU1;  IS  LANDING  COELETED 


CONTTNUl  DELAY:  'CSUTI'; 

ASSIGN:  A(2)  =  'CMTE'+26; 

ASSIGN:  A(3)  =  'CNLNE'+26+'UJEC'; 
ASSIGN:  'LTOE'  =  0; 

ASSIGN:  A(1)=TNCW; 

ASSIGN:  ’CLEN3'  =  C0(A(2)); 

ASSIGN:  ’SDIST’  =  C0(A(3)); 

ASSIGN:  ’CDIST’  =  'SDIST'; 

ASSIGN:  'TOES'  =  DP(20,1); 

ASSIGN:  'CTOE'  =  1:  NEXT  (CORRIDOR ) ; 


DELAY  CSUn 
INCREMENTS  A(2) 
INCREMENTS  A(3) 
ASSIGN  LTOE 
ASSIGN  TTO€ 

ASSIGN  CLEN3 
ASSIGN  SDIST 
ASSIGN  CDIST 
ASSIGN  TOES 
ASSIGN  CTOE  TTEN  GO 
TO  ELOCK  CCRRIDCR 


INITIALIZES  VALUES  WEEN  STARTING  AT  A  LEW  CUTHALL  DISTANCE  * 


LEWD  1ST  ASSIGN:  'LTOE'  =  TNCW  -  'TTOE' ;  ASSIGN  LTOE 

ASSIGN:  'CDIST'  =  'CDIST'  +  'INCRE';  ASSIGN  CDIST 

ASSIGN:  'TOES’  =  DP(20,1);  ASSIGN  TOES 

BRANCH,  1: 


IF,  'TUBES'  .Ed  0,  NEWDIST: 

ELSE,  C0NTTNU2; 

CCNTINL2  ASSIGN:  'CT1NE'  =  I; 

BRANCH,  1: 

IF,  'CDISP  .GT.  'CLENG',  NEWCCFR: 
ELSE,  CCRRIDCR; 


ANY  TURNS  AT  CURRENT 
CDIST 

ASSIGN  CXINE 


TESTS  IF  CCRRIDCR  IS 
CCNFLETED 


*  INITIALIZES  VALUES  WREN  STARTING  A  NEW  LAIRING  * 

jHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHt 


ICWlAfD 


ASSIGN:  'LNA6'  =  'LNUNB'  +  1; 
BRANCH,  1: 

IF,  'LRNC'.Ed  2,  LAND2: 
ELSE,  INISRED; 


ASSIGN  INLM3 


NEW  LANCING  CR  JCB 
CCNFLETED 


CCMfLETES  A  YARDING  CYCLE 

ASSIGN:  'CHINE'  =  'EETAO'  +  «BETA1'*'CDIST' 

;  ASSIGN  OT1NE 

DELAY: 

•CHINE'; 

DELAY  OTINE 

TALLY:  1, 

•CDIST' ; 

RECORD  CDIST 

TALLY:  2. 

'CHINE'; 

RECCRD  CHINE 

ASSIGN: 

'LOTIM'  =  WE  <  21 # 1 ) ; 

ASSIGN  LOTIM 

DELAY: 

'LOTIM' ; 

DELAY  LOTIM 

TALLY:  3, 

•LOTIM'; 

RECCED  LOTIM 

ASSIGN: 

•HOTIM'  =  WE(22,1); 

ASSIGN  HOTIM 

CELAY: 

'HOTIM'; 

DELAY  HOTIM 

T/LLY:  4, 

'HOilM'; 

RECCED  HOTIM 

ASSIGN: 

'LITTM'  =  CP(23,I); 

ASSIGN  LITIM 

DELAY: 

•LITIM'; 

DELAY  LITIM 

TALLY:  5, 

'LITIM'; 

RECORD  LITIM 

ASSIGN:  » 

ITINE'  =  'EETOO'  +  'BETlI'*'CDIsr» 

;  ASSIGN  ITINE 

DELAY: 

TINE' ; 

DELAY  ITTNE 

ASSIGN: 

'  HOTIM'  =  DP(24,1); 

ASSIGN  HTTIM 

DELAY: 

'  ID  i  1M' ; 

DELAY  IDTITJI 

TALLY:  6, 

'ITINE'; 

RECCFD  ITINE 

TALLY:  7, 

'HTTIM' ; 

RECCRD  HTTIM 

ASSIGN:  ' 

UHTIM'  =  RL(2S,I); 

ASSIGN  UHTIM 

DELAY:  ' 

UhTIM'; 

DELAY  UHTIM 

TALLY:  8, 

'UHTIM'; 

RECCRD  UHTIM 

ASSIGN:  'TCLEA'  =  DP(26,1); 


ASSIGN  TCLEA 


DELAY:  'TOEA'; 

TALLY:  9,  'TCLEA* ; 


DELAY  TOEA 
RECCFC  TOEA 


ASSIGN:  'CTM'  =  'CTM'  +  I; 
ERAfCH,  1: 

IF,  'CTM'  .LT.  'IMS',  CCFRIDCR: 
ELSE,  hECIST; 


ASSIGN  CTM 

TESTS  IF  ALL  TURNS 
CCNflFTED  AT  CDIST 


FINISH  GCLKT :  1, 1  :DISPCSE : 


ENDS  THE  SIMJLATICN 


CCUfN  1 
CCLINN  2 
OCLLHl  3 
CCIJ-N  4 

ccllnn  5 

CdlW  6 
OCLLW  7 
CCLINN  8 
CCLLNN  9 
CCLLWJ  10 
CCLIW  11 
COLON  12 


CrOE  NLM3ER 
OUTmj.  TINE  (SECONDS) 

OUmL  DISTANCE  (FEET) 

LATERE  CUimjL  TINE  (SECONDS) 
LATERAL  CUmi  TINE  (FEET) 
HOCKING  TINE  (SECONDS) 

LATERAL  INhWi  TINE  (SECONDS) 
Ihmi  TINE  (SECONDS) 

OKXXING  TINE  (SBXNDS) 

MNEER  CF  STEMS  LANCED 
TOTAL  CiOE  TINE  (SECONDS) 
COLON  WERE  DELAY  OCCURRED 


3 

23.00  100 

39.80 

30 

9.74 

6.20 

34.40 

16.09  1 

129.23 

4 

6.00  100 

57.80 

40 

14.85  240.35 

34.18 

1139  1 

364.77 

4 

5 

6.00  100 

39.00 

45 

14.28 

4.18 

32.75 

34.83  1 

131.04 

6 

10.00  120 

27.40 

50 

21.05 

45.17 

31.78 

15.29  1 

150.69 

8 

10.00  120 

27.70 

50 

32.82 

12.32 

32.63 

18.49  1 

133.96 

9 

11.00  120 

4330 

75 

13.67 

1639 

32.42 

17.12  1 

133.90 

10 

10.00  120 

29.20 

86 

24.99 

21.00 

32.78 

35.28  3 

153.25 

11 

11.00  120 

47.10  100 

26.02 

19.67 

32.27 

27.74  2 

163.80 

12 

11.00  120 

7230  100 

35.64 

44.75 

28.80 

16.79  1 

209.48 

13 

11.00  120 

26.50 

60 

26.14 

19.16 

32.55  124.60  1 

239.95 

6 

14 

12.00  150 

49.40 

20 

42.98 

16.08 

41.71  116.76  1 

278.93 

6 

15 

12.00  150 

35.40 

30 

81.72 

1631 

43.42  158.00  2 

346.75 

6 

16 

16.00  220 

16.92 

25 

48.99 

2838  163.00 

2035  1 

294.04 

5 

17 

18.00  220 

29.00 

35 

34.11 

15.17 

6333 

22.41  1 

182.22 

18 

16.00  220 

27.90 

35 

54.01  165.60 

57.68 

14.78  1 

335.97 

5 

19 

18.00  220 

29.00 

40 

24.44 

14.00 

5430 

25.09  1 

164.73 

20 

18.00  220 

34.80 

10 

31.62 

19.06 

60.32 

13.40  2 

177.20 

21 

16.00  220 

2830 

25 

82.07 

24.00 

61.46 

19.53  2 

231.56 

22 

17.00  220 

26.40 

30 

65.07 

2730 

61.19  233.37  1 

430.53  56 

23 

17.00  250 

17.80 

15 

40.03 

15.06 

78.42 

20.78  1 

189.11 

24 

21.00  250 

18.50 

20 

65.71 

15.85 

74.70 

38.42  2 

234.18 

26 

18.00  250 

29.30 

60  100.14 

26.18 

73.99  178.03  2 

425.64  46 

27 

22.00  325 

27.30 

20 

80.72 

10.82 

97.92 

32.91  1 

271.67 

3 

28 

22.00  325 

16.40 

25 

34.14 

13.73 

96. CB 

17.15  1 

199.45 

29 

22.00  325 

30.40 

30 

16.72 

13.43 

84.16 

27.57  1 

194.28 

30 

20.00  325 

16.30 

30 

25.85 

19.95 

82.45 

21.49  1 

186.04 

31 

21.00  325 

16.80 

25 

18.30 

17.19 

85.09 

28.81  1 

187.19 

32 

22.00  325 

42.10 

30 

29.67 

14.33 

91.77 

19.02  2 

218.89 

33 

22.00  325 

30.00 

50 

59.99 

16.63 

94.64 

173.65  1 

396.91 

6 

34 

21.00  325 

14.30 

60 

40.89 

16.03 

75.17 

10.92  1 

178.31 

35 

21.00  325 

36.80 

40 

20.01 

18.17 

79.25 

123.05  1 

298.28 

6 

36 

21.00  325 

27.60 

35 

32.74 

20.15 

78.66 

24.96  2 

205.11 

37 

20,00  325 

32.50 

75 

48.67 

20.00 

76.95 

32.53  2 

230.65 

38 

21.00  325 

28.00 

60 

72.78 

39.66 

80.63 

114.81  1 

356.88 

6 

39 

20.00  325 

30.90 

45 

31.87 

19.67  122.31 

19.14  1 

243.89 

40 

23.00  325 

52.60 

80 

35.06 

24.30 

89.74 

20.58  2 

245.28 

41 

21.00  325 

23.40 

50 

33.58 

16.93 

81.08 

28.53  2 

20432 

42 

21.00  325 

38.50 

60 

34.01 

29.45 

76.82 

26.99  2 

226.77 

43 

23.00  325 

28.30 

70  210.64 

18.03 

83.10 

30.10  3 

393.17 

3 

44 

21.00  325 

41J0 

70 

10.86 

32.38  144.86 

41.87  1 

292.27 

5 

45 

22.00  375 

24.20 

15 

33.67 

9.94 

97.10 

20.32  1 

207.23 

46 

23.00  375 

15.90 

10 

54.05 

10.67  102.21 

34.19  2 

240.02 

47 

25.00  375 

31.00 

50 

48.14 

10.78 

89.58 

15.15  1 

219.65 

46 

2S.00  375 

23.00 

45 

37.68 

12.66 

92.31 

15.24  1 

205.89 

49 

25.00  375 

17.80 

60 

48.13 

12.38 

91.98 

18.51  2 

213.80 

50 

2S.00  375 

48.80 

70 

29.3S 

30.48 

95.37 

18.15  1 

247.15 

52 

15.10  230  279.00 

10  350.00 

5.00 

59.51 

9.06  1 

717.67 

53 

18.10  250 

11.00 

2 

32.00 

61.00 

96.28 

14.74  1 

233.12 

54 

17.30  260 

22.00 

10 

7.00 

67.00 

96.38 

11.93  1 

221.61 

55 

17.50  290 

45.00 

1 

50.00 

68.00 

97.85 

9.39  1 

287.74 

4 

56 

17.60  280 

22.00 

4 

47.00 

74.00 

72.45 

19.49  1 

252.54 

50 

18.13  275 

15.00 

20 

25.00 

40.00 

74.03 

10. 13  1 

182.29 

59 

19.90  270 

13.00 

8 

28.00 

42.00 

73.24 

12.41  1 

18835 

60 

17.60  260 

16.00 

20 

51.00 

70.00 

79.71 

11.64  2 

245.95 

61 

18.50  260 

3.00 

25 

68.00  176.00 

68.98 

9.57  1 

344.05 

62 

17.20  260 

33.00 

35  115.00  176.00 

88.70 

23.25  2 

453.15 

63 

17.60  260 

21.00 

10 

45.00  163.00  178.55 

11.10  1 

436.2S 

5 

64 

20.00  260 

35.00 

30 

51.00 

73.00 

80.86 

6.87  1 

266.73 

65 

18.98  260 

23.00 

30 

44.00 

69.00 

76.41 

9.14  1 

240.53 

66 

20.10  260 

29.00 

35 

55.00 

80.00 

91.51 

7.97  1 

283.58 

67 

18.80  260 

45.00 

40 

83.00  109.00 

82.29 

1437  1 

352.46 

68 

18.50  260 

25.00 

20 

39.00 

68.00 

74.81 

8.87  1 

234.18 

69 

20.10  260 

25.00 

25 

64.00 

83.00 

73.00 

8.28  1 

273.38 

70 

18.40  260 

33.00 

35 

69.00 

92.00 

81.43 

11.37  1 

305.20 

71 

21.70  260 

28.00 

12 

55.00 

62.00 

78.95 

8.63  1 

254.28 

72 

18.50  265 

23.00 

20 

61.00 

77.00 

91.17 

20.64  1 

291.31 

73 

20.50  265 

46.00 

30  128.00  125.00  334.72 

16.15  2 

670.37 

5 

74 

17.20  265 

35.00 

10 

57.00 

76.00 

71.03 

9.54  2 

265.77 

75 

22.50  310 

39.00 

2 

44.00 

64.00 

82.80 

1132  1 

263.82 

76 

18.27  320 

38.00 

8 

65.00  134.00  107.83 

15.27  2 

378.37 

77 

21.96  320 

37.00 

20 

71.00 

95.00 

93.52 

36.62  1 

355.10 

78 

22.60  330 

5.00 

5 

24.00 

58.00 

97.50 

9.02  1 

216.12 

79 

19250  330 

10.00 

3 

38.00 

61.00 

96.94 

16.28  1 

241.72 

80 

23.30  330 

25.00 

20 

55.00 

82.00  282.12 

17.16  2 

484.58 

81 

22.30  330 

24.00 

25 

63.00 

90.00 

96.60 

9.61  1 

30531 

82 

21.90  340 

31.00 

30 

57.00 

85.00 

98.95 

24.69  1 

31834 

83 

22.50  340 

13.00 

5 

32.00 

52.00 

95.26 

11.24  1 

226.00 

84 

22.30  340 

27.00 

40 

53.00 

85.00  104.46 

10.02  1 

301.78 

85 

22.20  340 

35.00 

50 

54.00 

86.00  110.97 

31.16  1 

339.33 

66 

21.20  340 

30.00 

50 

81.00  110.00  166.03 

16.36  2 

424.59 

87 

23.96  340 

38.00 

45 

68.00 

99.00  107.79 

15.99  2 

352.74 

88 

20.10  340 

65.00 

50 

96.00  242.00 

92.54 

12.78  1 

528.42 

4 

89 

20.20  340 

42.00 

55  120.00  155.00  108.02 

16.23  2 

461.45 

90 

21.10  340 

67.00 

60  12S.00  159.00  105.75 

12.15  2 

490.00 

91 

21.20  345 

48.00 

65 

102.00  347.00  102.62 

15.58  2 

636.40 

4 

92 

19.40  345 

55.00 

70 

97.00  143.00  115.17 

12.50  2 

442. C7 

93 

20.30  345 

47.00 

10 

87.00  109.00  112.54 

10.74  1 

38638 

'v/ 

,  . 

vivIvVXvvOv 

s  s 


94  23.79  345  25.00  12  70.00  100.00  98.89 

95  22.31  345  37.00  25  80.00  107.00  103.27 

96  21.89  345  42.00  35  76.00  108.00  105.78 

97  21.64  345  46.00  70  58.00  130.00  119.63 

98  20.37  345  40.00  35  68.00  37.00  96.38 

99  20.78  345  34.00  70  75.00  105.00  103.91 

100  22.80  345  41.00  60  74.00  100.00  90.64 

104  13.00  300  48.00  70  45.00  17.00  40.00 

105  11.00  300  66.00  70  94.00  20.00  88.00 

106  13.00  300  77.00  70  73.00  890.00  159.00 

107  11.00  300  77.00  70  43.00  42.00  76.00 

108  20.00  300  79.00  60  14.00  13.00  36.00 

109  67.00  300  56.00  70  36.00  95.00  35.00 

110  36.00  300  66.00  70  62.00  32.00  65.00 

111  12.00  300  49.00  50  40.00  387.00  54.00 

112  12.00  300  41.00  40  34.00  518.00  133.00 

IB  73.00  300  47.00  50  112.00  84.00  42.00 

114  12.00  300  60.00  50  112.00  361.00  48.00 

115  B.OO  300  65.00  50  75.00  38.00  46.00 

116  97.00  300  68.00  70  108.00  371.00  55.00 

117  14.00  300  63.00  20  1Q3.00  21.00  44.00 

118  B.OO  300  49.00  20  107.00  10.00  66.00 

119  7.00  300  63.00  40  53.00  31.00  51.00 

120  164.00  300  76.00  40  19.00  31.00  44.00 

121  B.OO  300  62.00  50  62.00  31.00  56.00 

123  20.00  320  55.00  30  10.00  46.00  IB  .00 

124  15.00  350  49.00  50  54.00  25.00  35.00 

125  15.00  350  68.00  50  19.00  31.00  51.00 

126  14.00  350  88.00  50  10.00  25.00  42.00 

127  15.00  350  59.00  40  31.00  40.00  37.00 

128  16.00  350  59.00  40  37.00  38.00  46.00 

129  15.00  350  86.00  65  39.00  26.00  54.00 

130  16.00  350  66.00  50  35.00  129.00  62.00 

B1  17.00  350  56.00  60  31.00  16.00  78.00 

132  B.OO  75  17.00  10  27.00  5.00  37.00 

B3  B.OO  75  36.00  50  17.00  33.00  42.00 

134  10.00  160  23.00  20  48.00  10.00  33.00 

B5  15.00  160  19.00  10  20.00  15.00  43.00 

136  B.OO  160  18.00  20  20.00  184.00  91.00 

B7  47.00  160  80.00  30  22.00  23.00  43.00 

138  16.00  160  36.00  40  32.00  45.00  28.00 

139  16.00  160  23.00  50  25.00  53.00  68.00 

140  39.00  160  33.00  20  39.00  30,00  58.00 

141  11.00  210  87.00  20  52.00  86.00  191.00 

142  20.00  210  36.00  10  24.00  18.00  36.00 

143  21.00  210  69.00  4C  21.00  2S.OO  162.00 

144  20.00  210  67.00  40  31.00  33.00  107.00 

145  19.00  2 1C  28.00  30  17.00  25.00  87.00 

146  59.00  210  43.00  20  16.00  29.00  83.00 

147  22.00  210  38.00  20  30.00  39.00  65.00 

148  2C.00  210  25.00  30  61.00  12.00  66.00 


19.55  2 
B.73  1 
19.07  1 
11.70  1 
9.44  1 
8.15  1 
14.68  1 
25.00  1 
28.00  1 
142.00  1 
50.00  1 
24.00  1 
31.00  1 
17.00  1 
19.00  1 
15.00  1 
19.00  1 
17.00  1 
28.00  3 
22.00  3 
34.00  3 
20.00  1 
15.00  1 
17.00  1 
10.00  2 
19.00  1 
22.00  1 
25.00  1 
20.00  1 
28.00  1 
35.00  2 
20.00  1 
16.00  1 
18.00  1 
25.00  1 
43.00  1 
38.00  1 
17.00  1 
19.00  1 
19.00  1 
24.00  1 
11.00  1 
60.00  2 
21.00  1 
18.00  1 
89.00  2 
38.00  1 
61.00  1 
33.00  2 
43.00  1 
26.00  1 


337.23 

363.31 

372.74 

386.97 

273.19 

346.84 

343.12 

188.00 

307.00 

B54.00 

299.00 

186.00 

320.00 

278.00 

561.00 

753.00 

377.00 

610.00 

265.00 

721.00 

279.00 

265.00 

220.00 

351.00 

234.00 

263.00 

200.00 

209.00 

199.00 

210.00 

231.00 

240.00 

324.00 

216.00 

124.00 

184.00 
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